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ABSTRACT 


Force-limited vibration testing (FLVT) is effective in reducing the low-frequency vibra¬ 
tion test environment for CubeSats on the Naval Postgraduate School CubeSat Launcher 
(NPSCuL); however, the CubeSats are still subjected to high-frequency amplifications 
above 500 Hz from the NPSCuL structure. The excessive, high-frequency vibration has 
caused test failures and forces CubeSat developers to focus more on surviving environmen¬ 
tal testing instead of developing state-of-the-art technology. Whole-spacecraft isolation 
systems are often used to reduce these amplifications, but they currently exist only for large 
spacecraft and are too expensive to adapt for small satellites. These limitations motivated 
the combined use of FLVT and commercial-off-the-shelf (COTS) isolators on NPSCuL 
as a novel, practical, and low-cost method to reduce vibration levels for small satellites. 
This method significantly reduces the high-frequency amplification by up to 97%; the root- 
mean-square acceleration (Grms) over the entire test frequency range drops by up to 78%. 
These results should allow more sensitive and complex payloads to gain access to space 
on future NPSCuL missions and demonstrate how a worst-case environment on a small 
satellite can be improved. Implementing low-cost, COTS isolators on other small satellites 
and CubeSat launch applications could be useful as well. 
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CHAPTER 1: 

Introduction 


In recent years, small satellites have become increasingly sophisticated due to the collabo¬ 
ration and innovation that occur within the small satellite community. This is particularly 
true for CubeSats, which are 10-cm-square picosatellites of varying lengths; while it was 
more common for these satellites to be simple and robust when they were first built in 1999, 
more CubeSats now have state-of-the-art technology. As a result, the CubeSat components 
are more sensitive, and the satellites require less harsh launch environments than what was 
acceptable in the past. Primary spacecraft, which are typically large and expensive, are 
mounted at the top of the launch vehicle and have a controlled launch environment in the 
fairing. However, secondary payloads, such as CubeSats, are often mounted in unique lo¬ 
cations and orientations on the launch vehicle. The launch environment at these locations 
is not particularly well-controlled; the temperature range during flight is typically wider for 
a secondary payload compared to that of a primary payload, and the vibration environment 
is more harsh. One such mounting location for secondary payloads is on the aft end of the 
Centaur upper stage of the Atlas V launch vehicle. The Naval Postgraduate School CubeSat 
Launcher, or NPSCuL (pronounced NPS “cool”), carrying eight Poly-Picosatellite Orbital 
Deployers (P-PODs), is the most common secondary payload to fly in this location. 

1.1 Need for Reducing Vibration Levels 

As discussed in [1], the NPSCuL amplifies the response at the P-POD interface due to its 
cantilevered mounting configuration on the launch vehicle and thin-plate mounting walls, 
which results in a relatively high random vibration environment for the CubeSats flying on 
the Atlas V Aft Bulkhead Carrier (ABC) missions. While the ABC base input level of 7.6 
Grms [2] is considered reasonable for satellites, the Grms levels at the P-POD interface can 
be double that of the input. This amplification occurs even with force-limiting, which is a 
method of reducing the conservatism during a fixed-base test. These relative locations are 
highlighted in the vibration test configuration shown in Figure 1.1. 


1 





P-POD-to-NPSCuL 

Interface 


ABC Base Input 


Figure 1.1: ABC Input and P-POD Response Locations 


The resulting 15.1 Grms at the P-POD interface for a P-POD on NPSCuL is high compared 
to other P-POD launch interface environments, many of which are similar to NASA General 
Environmental Verification Standard (GEVS) vibration levels. In comparison, the GEVS 
levels are only 10.0 Grms [3] to satisfy the same test requirement. An acceleration spectral 
density (ASD) plot of a representative vibration profile on NPSCuL is shown in Figure 1.2, 
with the ABC input profile, indicating that there is a significant amount of energy across 
the entire test frequency range that needs be reduced. 
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Figure 1.2: ABC Input and P-POD Response Acceleration Profile Compar¬ 
ison 


Even though force-limited vibration testing has been used during ground testing for the past 
four missions, and is effective in lowering the low-frequency vibration test environment 
for CubeSats flying on NPSCuL, the high-frequency amplifications above 500 Hz were 
still problematic. The excessive, high-frequency vibration has caused test failures and has 
forced CubeSat developers to focus more on surviving environmental testing instead of 
optimizing satellite performance, thus motivating the need to reduce the vibration levels by 
re-designing NPSCuL, introducing isolators, or both [1]. 

1.1.1 Existing Spacecraft Isolation Systems 

Whole-spacecraft isolation systems are often used to reduce these amplifications, but they 
currently only exist for large spacecraft and are too expensive to adapt for small satellites. 
Moog CSA Engineering is the leading developer of whole-spacecraft isolation systems 
for large spacecraft, including the Softride and Shock Ring systems [4], These systems 
(shown in Figure 1.3) provide up to 50% attenuation of shock and vibration loads, but cost 
upward of $200k. Additionally, the SoftRide systems are not easily adaptable for payloads 
such as NPSCuL. They are designed to maintain the existing load path between the LV 
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and the spacecraft, so there is one isolator per interface fastener. The Softride Uniflex and 
Omniflex systems are also tuned to the fundamental frequency of the primary spacecraft, 
resulting in recurring costs for each mission due to the analysis required, in addition to the 
manufacturing costs. This increase in complexity during the spacecraft mate process to the 
launch vehicle, as well as the high cost, is not desirable for an auxiliary payload [5], [6]. 
The Shock Ring was adapted for the ESPA interface, but it was priced at $50k per ring, 
which is still relatively expensive for NPSCuL. Additionally, the attenuation was targeted 
for frequencies above 100 Hz, leaving any lower modes susceptible to amplification [7]. 



(a) SoftRide Uniflex (b) SoftRide Omniflex 



(c) SoftRide Shock Ring 


Figure 1.3: Whole-Spacecraft Isolation Systems. 

Source: “SoftRide Brochure.” [Online]. Available: 

http://www.csaengineering.com/literature/Space_Defense/Vibration Control/ 
MoogCSA_SoftRide_brochure.pdf 


Another limiting factor of adapting the Softride systems for small satellites is the surviv¬ 
ability temperature range of the viscoelastic material (VEM) used. VEM is a co m mon 
material used for vibration isolation; it has been incorporated in beam structures and other 
constrained-layer damping applications. [8], [9]. However, with a typical thermal operat¬ 
ing range of 10-35 °C (50-95 °F) and a survivability range of 5-50 °C (41-122 °F) [10], 
VEMs are considered sensitive to temperature for auxiliary payloads, for which the launch 
environment temperature can range from -45-38 °C (-50-100 °F). Therefore, using VEMs 
with constrained-layer damping was not considered for NPSCuL. 
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1.2 A Novel, Practical Method for Reducing Vibration 
Levels on Small Spacecraft 

Commercial-off-the-shelf (COTS) isolators are often used on launch vehicles and large 
spacecraft for avionics components and other sensitive equipment. However, these compo¬ 
nents usually weigh less than 22.7 kg (50 Ibm) and are mounted using four isolators, one 
at each corner of the component or mounting plate [11]. These COTS isolators have not 
been previously used for spacecraft similar to the size of NPSCuL (approximately 90.7 kg, 
or 200 Ibm) but are known to reduce high-frequency vibration. This dissertation discusses 
the combined use of FLVT and COTS isolators on NPSCuL and shows how a worst-case 
vibration environment can be reduced for small spacecraft in test and in flight. The result¬ 
ing performance is better than that of whole-spacecraft isolation systems designed for large 
spacecraft and is achieved at a small fraction of the cost. 

1.3 Background 

The development of the CubeSats, the P-POD, and NPSCuL is provided as background 
information. Unlike the larger satellites that have a direct interface with the launch vehicle, 
CubeSats require a deployer to serve as the launch vehicle interface as well as provide a 
means of deployment. Various CubeSat deployers exist, but P-PODs are the most common 
deployers used for the ABC missions. To maximize the CubeSat capacity for this secondary 
payload launch opportunity, NPSCuL is the launch vehicle interface between the ABC 
and up to eight P-PODs. The vibration test requirements for this system, which consists 
of CubeSats, P-PODs, and NPSCuL, are defined by the Test Requirements for Launch, 
Upper-Stage and Space Vehicles (SMC-S-016) [12]. The associated test levels to satisfy 
these requirements are also provided. 

1.3.1 The CubeSat Standard and P-PODs 

The CubeSat standard was developed in 1999 using available COTS components as a joint 
collaboration between California Polytechnic State University, San Luis Obispo (Cal Poly) 
and Stanford University. CubeSats are nanosatellites and picosatellites that are defined by 
a 10-cm cube weighing up to 1.33 kg; these are the specifications for a "lU" in the Cube¬ 
Sat Design Specification [13]. The most common other CubeSat sizes are L5U, 2U, and 
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3U, as shown in Figure 1.4. The original objective was to provide an educational platform 
for students to develop and launch satellites rapidly with low costs; the small size of these 
satellites allows for an entire mission life cycle to be completed within two years, or the 
average length of a master’s program [14]. Due to the small size and low cost of these satel¬ 
lites, CubeSats are considered secondary payloads. While universities are still provided the 
opportunity to participate in this unique educational experience, mainly through NASA’s 
Educational Launch of Nanosatellites (ELaNa) program, there is an increasing number of 
CubeSats that are developed by government organizations and private companies. 



1U SkHeton CAD Model RevD 



1.5U Skdeton 

CAD Model RevD 



2U Skeleton 
CAD Model 


RevD 



3U Skeleton 

CAD Model RevD 


Figure 1.4: Common CubeSat Sizes. 


Source: S. 

Naik, 

"CubeSat Constellation 

around 

Mars,” Dec. 

2013. 

[Online]. 

Available: 


http://ccar. colorado.edu/asen5050/proJects/projects_2013/Naik 
Siddhesh/Cubesats.html 


The P-POD was developed by Cal Poly and is the most commonly-used, standardized 
CubeSat deployment system in the world today. Shown in Figure 1.5, the P-POD is a 
non-separating, rectangular box made of aluminum, carrying any combination adding up 
to 3Us’ worth of CubeSats. The most common combinations are a single 3U or three lU 
CubeSats. The P-POD’s spring-loaded door is released by a non-explosive actuator (NEA), 
and the CubeSats are subsequently deployed by the main spring [13]. 
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Figure 1.5: P-POD and Main Spring. 

Source: CubeSat Design Specification, CubeSat Program, Cal Poly SLO 
Std., Feb. 2014. 


1.3.2 NPSCuL Overview 

The concept of clustering P-PODs in the NPSCuL was first published in 2007 in response 
to the increasing need for LF.S.-based launch opportunities [15]. The objective was to im¬ 
prove rideshare opportunities for the CubeSat community as a standard bus for interfacing 
multiple P-PODs to a launch vehicle. The standard bus was designed to be either a primary 
or secondary payload because late-stage manifesting was also considered in the event that 
a primary or secondary payload could not meet a milestone during the launch vehicle in¬ 
tegration flow [16]. The original NPSCuL was designed to maximize the space and mass 
available for an Evolved Expendable Launch Vehicle (EELV) Secondary Payload Adapter 
(ESPA) payload, which resulted in the D-advanced structure. This original design accom¬ 
modated 10 5U P-PODs, carrying up to 50 lU CubeSats [17]. 

In late 2008, NPSCuL was redesigned for the United Launch Alliance (ULA) ABC concept, 
which has a smaller envelope and mass allowable than the ESPA payload specification [18]. 
The ABC concept, shown in Figure 1.6, was conceived when ULA determined that one of 
the helium bottles on the aft end of the Centaur upper stage of the Atlas V vehicle was no 
longer necessary; this allowed for a secondary payload to take up the volume and mass. 
The ESPA 15-inch bolt-hole pattern is the standard interface on the ABC plate. 
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Centaur Aft 
Bulkhead 



Existing 

Mounting 

Brackets 


Aft Bulkhead 
Carrier (ABC) 


Auxiliary 
Payload (AP) 


Figure 1.6: Aft Bulkhead Carrier on Atlas V Centaur Aft Bulkhead. 
Source: Aft Bulkhead Carrier Auxiliary Payload User’s Guide, United 
Launch Alliance Std., May 2014. 


As shown in Figure 1.7, the resulting re-design of NPSCuL became the baseline NPSCuL 
design [19]. This design can accommodate eight P-PODs, four 6U dispensers, or various 
combinations of 3U and 6U dispensers. The 6U is a larger CubeSat form factor that is 
starting to be used; NPSCuL is compatible with both the dispensers built by Tyvak Nano- 
Satellite Systems, Inc. and Planetary Systems Corporation (PSC). The 6U dispenser built 
by Tyvak is derived from the Nanosatellite Launch Adapter System (NLAS) [20]. NPSCuL 
is an ESPA-compatible, five-sided aluminum structure with a non-separating adapter ring. 
The Splitter Auxiliary Device (SAD), an electrical interface box between the launch vehi¬ 
cle (LV) and the deployers, is also considered to be part of NPSCuL. Eight years after the 
NPSCuL concept was published, it has flown four times as part of the Operationally Unique 
Technologies Satellite (OUTSat), the Government Experimental Multi-Satellite (GEMSat), 
the Unique Lightweight Technology and Research Auxiliary Satellite (ULTRASat), and the 
Government Rideshare Advanced Concepts Experiment (GRACE) missions. These mis¬ 
sions were launched in 2012, 2013, and 2015, and a total of 46 CubeSats were successfully 
deployed. A fully-integrated payload is shown mounted to the ABC on the aft end of the 


8 








Atlas V Centaur stage in Figure 1.8. 


^AP 



Figure 1.7: NPSCuL Baseline Design 



Figure 1.8: OUTSat Mounted to ABC on Centaur Stage Aft End 
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1.3.2.1 UseofNPSCuL 

The NPSCuL geometry and interface to the LV is not unique for Evolved Expendeable 
Reusable Launch Vehicle (EELV) Secondary Payload Adapter (ESPA) payloads. A ring- 
to-box structure is common for many nanosatellites due to the ESPA standard, and the 
NPSCuL mass of approximately 90.7 kg (200 Ibm) is at the low end of the typical range 
of 90.7 to 227 kg (200 to 500 lbs) for an ESPA payload. FalconSat-3, mounted to one port 
of the ESPA ring, and the four ESPA payloads from the STP-S26 mission mounted to the 
Multiple Payload Adapter (MPA) at the forward end of the Minotaur IV launch vehicle are 
shown in Figure 1.9. These payloads did not have any reported vibration test issues, and the 
STP-S26 test environments were much less harsh due to their primary payload status, but 
they are examples of other ESPA payloads with similar geometry and boundary conditions 
between the launch vehicle and components mounted to the primary spacecraft structure. 



Figure 1.9: FalconSat-3 (left) and ESPA Payloads on STP-S26 (right). 
Sources: “Falconsat-3," 2002. [Online]. Available: 

https://direct 0 ry.e 0 p 0 rtal. 0 rg/web/e 0 p 0 rtal/satellite-missi 0 ns/f/ 
falconsat-3; “Stp-s26,” 2002. [Online]. Available: 

https://directory, eoportal.org/web/eoportal/satellite-missions/s/stp- 
s26 
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1.3.3 Definition of Test Levels 

Satellite vibration test requirements can come from various sources and often depends on 
the responsible organization. As a a result, the definition of test levels can also vary. The 
maximum predicted environment (MPE) at the ABC is the base input shown in Figure 1.2; 
the MPE is also considered the acceptance level for vibration testing. Following the guide¬ 
lines in [12], protoqualification and qualification levels, defined as MPE -i-3 dB and MPE +6 
dB, respectively, are used for different reasons throughout the spacecraft development cy¬ 
cle. Typically, the qualification levels are used on an engineering development unit (EDU) 
to show that the design has adequate margins for surviving launch loads. The acceptance 
levels are used for testing performed on the flight unit in its final configuration after a 
flight-like EDU has been tested to qualification levels; these levels are meant to show that 
the unit has been properly built, which also demonstrates process control for multiple flight 
units. The protoqualification levels, which are performed on a flight unit, are used on a 
case-by-case basis when there is no EDU available or when the risk of testing an EDU to 
qual is higher than testing a flight unit with extra margin to gain confidence in its structural 
integrity. With this approach, the flight unit must also be subjected to an acceptance test. 
In both the qualification and protoqualification test approaches, subsequent flight units can 
then be tested to acceptance levels [12]. It is noted that the protoqualification approach is 
not the same as the protoflight approach, even though the terms are sometimes used in¬ 
terchangeably. In the protoflight testing strategy, the flight unit is only subjected to the 
protoqualification test levels, and it is not re-tested at the acceptance level. This approach 
is often used to reduce hardware cost, but it carries schedule and technical risk [21]. 

1.4 Dissertation Summary 

The progression of methods used to reduce the vibration environment on NPSCuL, which 
include the benefits and shortcomings of using FLVT and the isgorid design, are presented 
in Chapter 2. In Chapter 3, the driving design requirements, the isolator selection, and 
the associated re-design of NPSCuL to accommodate each isolator are discussed along 
with a summary of the test configurations, test inputs, and test set-up. The test results 
are presented and discussed in Chapter 4, as well as the performance assessment of both 
isolators. Chapter 5 describes the finite element models that were constructed and used 
for analysis for each test configuration, as well as the loads analyses that were performed 
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on the configuration that is the most likely to fly on a future ABC mission. Concluding 
remarks, the summary of contributions, and future work are provided in Chapter 6. 
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CHAPTER 2: 

Evolution of Vibration Reduction Methods on 

NPSCuL 


Changing the test set-up, the structure itself, and the test levels from the LV provider were 
considered as methods to reduce the vibration levels on NPSCuL. Only changes to the test 
set-up and the NPSCuL design will be discussed; at this time, LFLA is performing a study 
to reduce the test levels, but updated test levels are not available. [1]. 

2.1 Force-Limited Vibration Testing 

The qualification and acceptance random vibration tests performed at the system level for 
the four NPSCuL missions to date all used force-limited vibration testing (FLVT) [19]. It 
is based on the concept that due to the vibration absorber effect, there is a physical limit 
on the force at the base of the test article that is often exceeded at the test article’s resonant 
frequencies with a fixed-base test set-up. Similar to the acceleration input, the force limit 
represents the envelope of the total expected force at an interface during flight [22]. The 
resonances at the NPSCuL fundamental frequencies are effectively reduced with the force- 
limited test set-up, resulting in a notched acceleration input that is less conservative than the 
prescribed acceleration input in [2] but still satisfies the ULA vibration test requirements. 
Traditionally, the test article is mounted directly to the shaker adapter plate for a fixed-base, 
random-vibration test. The input is only controlled with respect to acceleration. This type 
of single-channel control is shown in Figure 2.1. 
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^ Acceleration Input (Reference) 


Acceleration Input 
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Acceleration Response \ 


Acceleration, Measured \ 




Acceleration Response 


Test Article 


j error Acceleration, Measured N- 


error = (Acceleration Input) - (Acceleration, Measured) 


Figure 2.1: Single-Channel Control 


In an FLVT set-up, such as the one used for NPSCuL shown in Figure 2.2, dual control 
is used; both the acceleration input and the force limit are the control channels. The force 
measured between the shaker and the base of the test article is compared to the force limit 
while the acceleration input is compared to the reference input; both are done simultane¬ 
ously by the shaker controller software. The software nominally controls the shaker using 
the acceleration input, but when the force limit is met or exceeded, the control switches to 
the force limit, which is set as a watchdog. At these frequencies where the force is being 
controlled, the resulting acceleration is only being measured. A top-level depiction of the 
dual-control loop used in an FLVT set-up is shown in Figure 2.3. 
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Figure 2.2: FLVT Set-Up of NPSCuL, X-Axis 
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Figure 2.3: Dual-Control Loop 


2.1.1 FLVT Background 

FLVT was developed at the Jet Propulsion Laboratory (JPL) for large spacecraft and Space 
Shuttle payloads as a way to reduce over-test resulting from a traditional, acceleration- 
controlled, fixed-base vibration input [23]. It is well known that the rigid fixed-base bound¬ 
ary condition is conservative and not representative of the flight condition, but the margin 
associated with vibration testing is also widely accepted by the aerospace industry. This 
type of set-up can be represented as shown in Figure 2.4. 
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Base Excitation (Prescribed Motion) 

Unnotched Acceleration, Ao X2t^2t ^2 



[Moving Support (Vibration Table) 


Figure 2.4: Fixed Base Test Set-Up Represented as 1 DOF Flarmonic Os¬ 
cillator Model with Prescribed Support Motion. 

Adapted from NASA-HDBK-7004B - Force Limited Vibration Testing, 
NASA Std., Jan. 2003. 


In flight, however, the spacecraft is mounted to a flexible structure, not a rigidly fixed struc¬ 
ture; therefore, the system can be simplified to a two degree-of-freedom (2DOF) harmonic 
oscillator model with prescribed motion, as shown in Figure 2.5. 


Base Excitation (Prescribed Motion) 

Urtnotched Acceleration, Ao %, J 

_c: 

flC Plate 
(Source) 


Interface Force 
I and Accleration ) 






2x ^2f ^2 


spacecraft 

(load) 

m2 




Moving Support (Launch Vehicle) 






Figure 2.5: FLVT Set-Up Represented as 2DOF Harmonic Oscillator Model 
with Prescribed Support Motion. 

Adapted from NASA-HDBK-7004B - Force Limited Vibration Testing, 
NASA Std., Jan. 2003. 
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This set-up results in the vibration absorber effect, in which the second mass, or the load 
(m 2 ), acts as a tuned absorber such that the anti-resonance frequency of the first mass, 
or the source (mi), is at the fundamental frequency of the load. The anti-resonance is a 
frequency at which the response of the excited mass is zero for an undamped system, or 
small for a damped system. The drop in magnitude of the source acceleration, at the fixed- 
base fundamental frequency of the load, is the theory behind the resulting notch in the 
acceleration control input when using a force limit. The transmissibility equations used to 
show this effect using NPSCuL in Section 2.1.1.3 are derived as follows. 

2.1.1.1 Transmissibility Equation - Forced Excitation of SDOF System Due to Sup¬ 
port Motion 

The transmissibility equation for the SDOF mass due support motion is derived from the 
equation of motion for a mass-spring-damper system with prescribed support motion using 
absolute coordinates, where x is the coordinate of the lumped mass, ^ is the damping ratio, 
and (On is the natural frequency of the system. 

i:' -I- 2^ (OnX + (OnX = + 2l^(Ony (2.1) 

The prescribed motion, y, is a function of time, where Q is the excitation frequency. 

y{t) = Ysin{Q.t) (2.2) 

To solve the ordinary differential equation, the particular solution for x(t) and y(t) are as¬ 
sumed to be as follows, where X and Y are the amplitudes of displacement of the lumped 
mass and the support motion, respectively. 

x{t)=Xe^^^ (2.3a) 

y{t) = YeJ^ (2.3b) 
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Substituting Equation 2.3 into Equation 2.1 and solving for the ratio between X and Y 
yields the non-dimensionalized transmissibility equation shown in Equation 2.4. 



2.1.1.2 Transmissibility Equation - Forced Excitation of 2DOF System Due to Sup¬ 
port Motion 

The transmissibility equations for each of the masses with prescribed base motion can be 
derived starting with the equations of motion for the 2DOF system in matrix form are 
shown in Equation 2.5, assuming ci is zero [24]). 

m 0 x\ ^ C2 -C2 Xi ^1+^2 -h Xi _ coly 

0 m2 X2 -C2 C2 X2 -^2 ^2 X2 0 

Similar to the SDOE case, the prescribed motion is defined in Equation 2.2. The frequency 
response functions (ERE) between the base and each mass. Hi and H 2 , can be calculated 
by taking the inverse of the impedance matrix, which is defined in Equation 2.6. 

[Z(Q)] = [K- Q?M + ;QC] (2.6) 

The resulting FRFs, calculated using Matlab (see Appendix D), are shown in Equation 2.7. 


Hi 

H2 


j{m2^ - k2) + cQ. 

j{—kik2 + kim2ki^ + k2miQ?- + k2m2kl^ — mim2kl^) + {c2kiQ.-\- C2miQ? -\-c 2 m 2 O?) 

(2.7a) 

_ k2 + jC2^ _ 

{kik2 — kim2^ — k2niikl^ — k2m2Q?'-\-mim2Q.‘^) + j{c2kiQ, —C2miQ.^ —C2m2Q.^) 

(2.7b) 


19 



Equation 2.7 is now in the form [25] 


Hi = 


A + jB 
C + jD 


and can be transformed into the form shown in Equation 2.9. 


Hi 


A2+52 

C2+E)2 


( 2 . 8 ) 


(2.9) 


Non-dimensionalizing Equation 2.7 in the form of Equation 2.9 yields 


Ai 



_ (2Cri^)2 + (r2-g2)2 _ 

(1 + {\+ll)q^)rl + q^Y + (2Cri^)2(l - r2(l +^))2 


A2 = 



(1 + {\+ll)q^)rl + q^Y + (2Cri^)2(l 


r 2 ( l +^))2 


(2.10a) 

(2.10b) 


using the ratios shown in Equation 2.11 [24]. 


a 


r\ = — 

(2.11 a) 

©1 


Q. 

r2 = — 

(2.IIb) 

(02 


m2 

B = — 

(2.11 c) 

m\ 


(02 

q = — 

(2.1 Id) 

©1 



Equations 2.10a and 2.10b are the transmissibility equations for the 2DOE system, similar 
to Equation 2.4 for the SDOE system. 

2.1.1.3 Force Limit Example with NPSCuL as a Single Degree-of-Freedom 

If NPSCuL is assumed to be a simple lumped mass weighing 85.7 kg (189 Ibm), the fixed 
base set-up can be modeled using the mass of NPSCuL and an assumed stiffness. In this 


20 



example, the stiffness is chosen such that the fundamental frequency of NPSCuL in the 
SDOF model is 117 Hz, as shown in Figure 2.6. This fundamental frequency corresponds 
to the first mode of NPSCuL in the axial (Z) direction. When the test set-up is modeled as 
a 2DOF system, the anti-resonance of the ABC plate is at 117 Hz, which is a result of the 
vibration absorber effect. 



Figure 2.6: Transmissibility of NPSCuL (SDOF Model) and ABC Plate 
(2DOF Model) 


The more realistic acceleration input corresponding to the acceleration of the ABC plate 
(mi), near the fundamental frequency of NPSCuL, is then predicted using the transmissibil¬ 
ity of the ABC plate in the 2DOF model. The mass of the ABC plate is estimated to be 9.07 
kg (20 Ibm). The transmissibility of the ABC plate is then multiplied by the un-notched 
acceleration input (corresponding to the support motion, A^, in Figure 2.5) to obtain the 
predicted notched acceleration input. In this example, the random vibration environment 
defined at the ABC interface, shown in Figure 1.2, is used as the un-notched acceleration 
input. 
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As shown in Figure 2.7, the actual acceleration input used during test compares well with 
the notch predicted in the 2DOF model. The two peaks in the predicted acceleration in¬ 
put are irrelevant; force limiting is only intended to reduce the input at the fundamental 
frequency of the system, so the 2DOF model is only used to validate the resulting notch 
in the FLVT test set-up. Additionally, the force limit has been shown to be conservative 
compared to actual flight data [23]. 



Figure 2.7: Predicted and Actual Force-Limited Acceleration Input 


2.1.2 Force Limit Derivation Methods 

The force limit can be determined using one of two analytical derivation methods, the 
simple and complex two-degree-of-system (TDFS) methods, or the semi-empirical method. 
The analytical methods require more knowledge of the system and have been shown to be 
more conservative than the semi-empirical method. All three methods will be described, but 
due to the unavailability of residual and modal masses of the launch vehicle interface and 
the trend that the analytical methods are typically more conservative, the semi-empirical 
method was used for NPSCuL. An example of how the force limit is calculated using the 
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simple and complex TDFS methods can be found in [23] and is demonstrated for NPSCuL 
here for all three methods. 

2.1.2.1 Simple Two Degree-of-Freedom System Method 

The simple TDFS method of deriving force limits is also referred to as the frequency shift 
method of predicting maximum force [26]. In this method, the following equation for 
a single-output, single-input force spectral density, is likened to F=ma, where Sff is the 
force spectral density, M 2 is the load apparent mass, and Saa is the acceleration spectral 
density [27]]. 

S//(Q) = |M2(n)|2Sa«(n) (2.12) 

The apparent mass is an FRF between the force and acceleration, and all three quantities 
are functions of frequency. The force and acceleration spectral densities are defined at the 
interface between the source and the load (see Figure 2.5). The dynamic amplification 
equation for the force spectral density, otherwise known as the normalized force limit, is 
simply the square of the right-hand side of Equation 2.4, where ^ is replaced with Q = ^ 
[27]. 

%(n) i + iff 

Q is also known as the quality factor and corresponds to the transmissibility at resonance 
for an SDOF system [21]. Note that M 2 is the scalar load mass, which is not the same 
as the apparent mass in Equation 2.12. For a classical undamped absorber, the squared 
ratio between the excitation frequency and the natural frequency of the system is shown in 
Equation 2.14, where /r is the mass ratio, ^ [25]. 

(s) + (2.14) 

Equation 2.14 is substituted into Equation2.13 to generate Figure 2.8. This plot is used to 
determine the maximum normalized force spectral density for varying mass ratios. These 
normalized values also correspond to the value in the semi-empirical method. 
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Figure 2.8: Normalized Force Spectral Limit from Simple TDFS Method. 
Source: NASA-HDBK-7004C - Force Limited Vibration Testing, NASA 
Std., Nov. 2012. 


The actual force limit is determined by multiplying the normalized value by the acceler¬ 
ation spectral density value; this is done at frequencies throughout the force-limited fre¬ 
quency range. It is recommended to determine the ratio of load to source masses using the 
asymptotic apparent masses [28]. It is also recommended that the force-limited frequency 
range be split into one-third octave bands, which can be calculated as shown in Equation 
2.15 [29]. 


fn+l 


fn 


= 2 * 


(2.15) 


Using the masses of NPSCuL and the ABC plate, the mass ratio is approximately 10. 
Therefore, the maximum normalized force spectral density is approximately 1.4 from Fig¬ 
ure 2.8. The one-third frequency bands, the ASD values from the ABC levels at MPE 
corresponding to the center frequencies, and the resulting force limit values are shown in 
Table 2.1. 
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Table 2.1: Force Limit Calculated Using Simple TDFS Method for NPSCuL 


Frequency Range Split Into 1/3 Octave Bands 

Acceleration 

Force 

Low 

Center 

High 

Spectral 

Spectral 

Frequency, 

Frequency, 

Frequency, 

Density, 

Density, 

Hz 

Hz 

Hz 

Saa ig^/Hz) 

Sff (N^Hz) [Ibf/Hz] 

19.8 

25 

31.5 

0.0515 

1.16E+04 [1.52E + 03] 

31.7 

40 

50.4 

0.125 

2.81E+04 [2.60E + 03] 

50.8 

64 

80.6 

0.125 

2.48E+04 [6.3 IE+ 03] 

81.0 

102 

128.5 

0.125 

1.31E+04 [5.57E + 03] 

128.6 

162 

204.1 

0.125 

1.31E+04 [2.95E + 03] 

204.8 

258 

325.1 

0.1101 

1.15E+04 [2.95E + 03] 

325.4 

410 

516.6 

0.0487 

5.10E+03 [2.59E + 03] 

396.9 

500 

630.0 

0.003 

3.14E+02 [1.15E + 03] 


2.1.2.2 Complex Two Degree-of-Freedom System Method 

The expansion of the simple TDFS model to account for the modal (m„) and residual (M„) 
masses, for both the source (ABC plate, m\ and M\) and the load (NPSCuL, m 2 and M 2 ), 
is shown in Figure 2.9. The modal mass refers to the effective mass of a fundamental 
mode, and the residual mass refers to the sum of the effective masses that are above the 
fundamental mode; these values are usually obtained from a finite element model (FEM). 
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Figure 2.9: Complex TDFS Method Model. 

Source: Y. Soucy, V. Dharanipathi, and R. Sedaghati, “Compar¬ 
ison of methods for force limited vibration testing," in Proceed¬ 
ings of the IMAC XXIII Conference, vol. 31. 31 January3 

February 2005, Orlando Fla. Society for Experimental Mechanics, 
Bethel, Connecticut, Paper No. 295, 2005. [Online]. Avail¬ 
able: http://sem-proceedings.eom/23i/sem.org-IMAC-XXIII-Conf-s03p02- 
Application-Operating-Data-Scaling-Techniques-Using-Multiple.pdf 


In the complex TDFS method, the modal and residual masses are accounted for in the 
H2{0.) term of Equation 2.12, as shown in Equation 2.16. 


\Tj ro'i |2 “^2)+ ^2)^+ ^2(1 /Ql 

= mwTwm 

The variables in Equation 2.16 are defined in Equation 2.17. 

P^ = -^B±^Vb^-4C 

(I + B 4-CCl)/{o>2/(OlY + (I + H + ^ 1 (^ 2 ) 

ts = - 


c = 


1 + /i 

\ ^ ji ^ ulphd 1 -|- jU (X2 


(l+ju)(a)2/a)i)2 
a\ = m\/M\ 

0 C 2 = m2/M2 

II = M2/M1 


(2.16) 


(2.17) 

(2.18) 

(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 
( 2 . 22 ) 


However, the peak acceleration and peak force cannot be assumed to occur at the same 
frequency; therefore, the acceleration peak ratio between ai to the external acceleration of 
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the residual mass, Ai, must be calculated to appropriately scale H 2 { 0 .), where Hi (Q) is the 
FRF between the source and the drive point. [26]: 


a\ 

2 

Hi 

Ax 


H1+H2 


(2.23) 


Both the load and source are accounted for in this model, so the ratio of the resonance 
frequencies between them is also a variable in this force limit derivation method. Therefore, 
the normalized maximum force limit values are tuned to ensure that the maximum value 
was found for all mass and damping combinations considered [26]. Due to the greater 
number of ratios involved in the complex TDFS method, the normalized force spectral 
limit is tabulated in terms of tti, a 2 , and fx. These values, for a Q of 20, are shown in 
Table 2.2 [28]. 
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Table 2.2: Normalized Force Limit Spectrum for Complex TDFS with Q = 

20 . 

Source: NASA-HDBK-7004C - Force Limited Vibration Testing, NASA 
Std., Nov. 2012. 
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When calculating the force limit using the complex TDFS method for NPSCuL, the modal 
and residual masses of NPSCuL were obtained from the FEM. The modal and residual 
masses of the ABC plate are unknown, so these values were estimated. Theses masses, the 
resulting ratios, the corresponding normalized force spectral density from Table 2.2, and the 
actual force spectral density for each one-third octave band are summarized in Table 2.3. 
The same octave bands and corresponding ASD values used in the simple TDFS method 
are used for this calculation as well; therefore, only the center frequencies of each band are 
shown here for reference. 
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Table 2.3: Force Limit Calculated Using The Complex TDFS Method for 
NPSCuL 


Center 

Freq.,Hz 

Modal 

Masses, 
kg [Ibm] 

Residual 

Masses, 
kg [Ibm] 

Mass 

Ratios 

Normalized 

Force 

Spectral 

Density 

Force 

Spectral 

Density 

(%)’ 

N^/Hz 

[IbfiHz] 

m\ 

m2 

Ml 

M2 

a\ 

052 

20 

0.00 

[0.00] 

0.00 

[0.00] 

9.08 

[20.0] 

85.8 

[189] 

0.00 

0.00 

1 

4.77E+03 

[1.52F + 03] 

25 

0.00 

[0.00] 

0.00 

[0.00] 

9.08 

[20.0] 

85.8 

[189] 

0.00 

0.00 

1 

8.18E+03 

[2.60F + 03] 

40 

2.27 

[5.00] 

0.03 

[0.07] 

6.81 

[15.0] 

85.78 

[189] 

0.25 

0.00 

1 

1.98E+04 
[6.3 IF+ 03] 

64 

6.72 

[14.8] 

5.18 

[11.4] 

0.09 

[0.20] 

80.60 

[178] 

0.74 

0.06 

4 

7.02E+04 

[5.57F + 03] 

102 

0.00 

[0.00] 

22.00 

[48.5] 

0.09 

[0.20] 

58.6 

[129] 

0.00 

0.26 

6 

6.59E+04 
[2.95F + 03] 

162 

0.00 

[0.00] 

0.00 

[0.01] 

0.09 

[0.20] 

58.6 

[129] 

0.00 

0.00 

1 

1.99E+04 

[2.95F + 03] 

258 

0.00 

[0.00] 

0.00 

[0.00] 

0.09 

[0.20] 

58.6 

[129] 

0.00 

0.00 

1 

1.75E+04 

[2.59F + 03] 

410 

0.00 

[0.00] 

0.04 

[0.09] 

0.09 

[0.20] 

58.6 

[129] 

0.00 

0.00 

1 

7.73E+03 
[1.15F + 03] 

500 

0.00 

[0.00] 

0.00 

[0.00] 

0.09 

[0.20] 

58.6 

[129] 

0.00 

0.00 

1 

4.77E+02 

[7.06F + 01] 


It is noted that if the analytical approach is used, the force limit should be calculated using 
both the simple and complex TDFS methods, and the force limit should consist of the max 
values [23]. Also, modal and residual mass values of the launch vehicle are not typically 
available. 
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2.1.2.3 Semi-Empirical Method 

The semi-empirical method is based on the extrapolation of interface force data measured 
during flight on various launch vehicles and mounting structures. The random vibration 
equations used in the semi-empirical method to calculate the force spectral limit are similar 
to 2.12 and are shown in Equation 2.24, where S/f is the force spectral density, Q is the 
driving frequency, M 2 is the total mass of the load, Saa is the acceleration spectral density, 
tao is the fundamental frequency of the load, and n is a positive constant. 

%(a) = C^MlSaai^) a < too (2.24) 

%(a) = c 2 M| 5 a«(a)/(too/a)” a > omeguQ (2.25) 

n is usually 2, but can be adjusted to fit experimental measurements of the payload apparent 
mass, which is the FRF between the measured force and acceleration. The term can be 
obtained from Figure 2.8, but using a value of 2 has also been shown to envelope the 
interface force data. Additionally, C is usually no higher than 1.4 due to the vibration 
absorber effect [27]. The decrease in payload residual mass as a function of frequency is 
accounted for the in the ///o term [23], [28]. 

The force limit calculated using the semi-empirical method for NPSCuL is shown in Fig¬ 
ure 2.10, along with the the force limits calculated using the simple and complex TDFS 
methods. For the semi-empirical method, M 2 is the mass of NPSCuL, and too is 117 Hz. 
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Figure 2.10: Force Limit Calculations for NPSCuL 


Consistent with the force limits calculated for various spacecraft and spacecraft compo¬ 
nents in [23], [28], the force limit calculated using the analytical methods (max values of 
both simple and complex TDFS methods) is generally more conservative than the semi- 
empirical method. 

2.1.3 Implementation of Force Limiting on NPSCuL 

Examples of a force limit for NPSCuL, derived using the semi-empirical method, and the 
resulting acceleration input at the base of NPSCuL are shown in Figure 2.11 and Fig¬ 
ure 2.12. In Figure 2.11, the measured force is expected to hit the force limit near the 
fundamental frequencies; it is not expected to hit the limit at other frequencies because of 
the absence of resonances. Because the modal and residual masses of the source are not 
typically readily available, the semi-empirical method is used to derive the force limit [23]. 
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Figure 2.11: Example Force Limit to 500 Hz, Z-Axis 


The notches in the acceleration input, shown in Figure 2.12, are near the frequencies at 
which the measured force does hit the force limit because the acceleration input has been 
adjusted according to the differences in the measured forces and the force limit at those 
frequencies. At all other frequencies (i.e., where the measured force does not meet the 
force limit) it is expected that the force-limited input is the same as the unnotched input. 
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Force Limited Input ^^«Unnotched Input Fundamental Frequency 


Figure 2.12: Resulting Notch at Fundamental Mode in Acceleration Input, 
Z-Axis 


The X-axis response at the P-POD #2 interface shown in Figure 2.13 and the corresponding 
Grms values shown in Table 2.4 highlight both the effectiveness of force-limiting and the 
inadequacies that still exist without further modifications to NPSCuL. The overall Grms is 
reduced by 56%, but the responses, shown in the frequency domain, are a representative 
case of the amount of high-frequency content that is present, even if an FLVT set-up is 
used. 
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—Unnotched Input FLVT Input—P-POD 2 Response, Z-Axis 


Figure 2.13: X-Axis Response at P-POD ^2 Interface (P2M2 Data at 
MPE +0 dB) 


Table 2.4: P-POD 2 Response, with and without FLVT 


Level 

Base Input Grms 

P-POD Interface Grms 

With Force Limiting 

6.18 

11.8 

Without Force Limiting 

7.6 

26.7 


It is noted that there is no relief past the force-limited range; however, test data shows that 
there is a significant amount of energy above 500 Hz at the P-POD interface, which is the 
frequency range of concern to many CubeSats. However, force limiting is only intended to 
reduce the input at the fundamental frequency of the system; it is not meant to change the 
properties of the system. 
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2.2 Isogrid NPSCuL 

FLVT is intended to change the test set-up by reducing the acceleration input at the fun¬ 
damental frequency, and it is very effective for that purpose on NPSCuL. However, FLVT 
is not intended to be a mechanism for changing the resonant properties of a system, so it 
was recognized that NPSCuL had to be modified to reduce the vibration environment at the 
P-POD interface. The main intent of the design changes were to stiffen the structure so that 
the displacement, and energy required to excite the resonant modes, would decrease. For a 
system subjected to a static force, the displacement decreases when the stiffness increases; 
in a dynamic environment, more energy is required to excite higher modes. The baseplate 
and wall thicknesses were increased, and an isogrid design was implemented to keep the 
primary structure mass to a minimum. Additionally, the baseplate and adapter ring were 
combined into a single part, called the unibase. The resulting re-designed isogrid NPSCuL 
structure is shown in Figure 2.14 and Figure 2.15 [19]. 



Figure 2.14: Isogrid NPSCuL Design. 

Source: V. Kaushish, “Force limited vibration testing and subsequent re¬ 
design of the naval postgraduate school CubeSat launcher,” Master’s thesis, 
Monterey, California: Naval Postgraduate School, 2014. [Online]. Avail¬ 
able: http://calhoun.nps.edu/handle/10945/42656. 
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Figure 2.15: Unibase. 

Source: V. Kaushish, “Force limited vibration testing and subsequent re¬ 
design of the naval postgraduate school CubeSat launcher,” Master’s thesis, 
Monterey, California: Naval Postgraduate School, 2014. [Online]. Avail¬ 
able: http://calhoun.nps.edu/handle/10945/42656. 


The isogrid NPSCuL design was built and tested, and test results were compared with the 
baseline NPSCuL vibration test results. It was found that, although the first fundamental 
frequencies of the isogrid NPSCuL design are higher than the baseline NPSCuL fundamen¬ 
tal frequencies, FLVT is less effective on the isogrid NPSCuL structure. This is because 
the force roll-off frequency associated with FLVT, which coincides with the fundamental 
frequency of the system, begins at a lower frequency when testing the baseline, less-stiff 
NPSCuL, allowing for a broader roll-off frequency range and slightly deeper notches. This 
is depicted in Figure 2.16, where the notch at 45 Hz, the fundamental frequency of the base¬ 
line NPSCuL, is deeper than the notch around 60 Hz, which coincides with the fundamental 
frequency of the isogrid NPSCuL in the X-axis. 
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Figure 2.16: X-Axis Notched Acceleration Control, Baseline and Isogrid 
NPSCuL at MPE -3 dB. 

Source: V. Kaushish, “Force limited vibration testing and subsequent re¬ 
design of the naval postgraduate school CubeSat launcher,” Master’s thesis, 
Monterey, California: Naval Postgraduate School, 2014. [Online], Avail¬ 
able: http://calhoun.nps.edu/handle/10945/42656. 


The resulting overall Grms levels measured at the P-POD-to-NPSCuL interface are higher 
for the stiffer isogrid structure, even when utilizing FLVT. This is shown in Figure 2.17, 
which contains the X-axis data taken at P-POD #2 for both the baseline and isogrid 
NPSCuL designs. 
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Figure 2.17: X-Axis Random Vibration Response, Baseline and Isogrid 
NPSCuL at MPE -3 dB. 

Source: V. Kaushish, “Force limited vibration testing and subsequent re¬ 
design of the naval postgraduate school CubeSat launcher,” Master’s thesis, 
Monterey, California: Naval Postgraduate School, 2014. [Online]. Avail¬ 
able: http://calhoun.nps.edu/handle/10945/42656. 


The results of the isogrid NPSCuL testing show that stiffening a structure may not reduce 
vibration test levels when utilizing FLVT. Additionally, there is no relief in amplitude at 
the P-POD interface, especially above 500 Hz. It was determined that either introduction 
of damping or utilization of vibration isolation techniques would probably be required to 
achieve better vibration environments [1], [19]. 
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CHAPTER 3: 

Incorporating Isolators on NPSCuL 


Stiffening the structure did not produce favorable results even when combined with FLVT, 
so incorporating isolators on NPSCuL was investigated as a way of reducing the stiffness to 
improve the vibration environment. The lower stiffness would potentially result in a lower 
roll-off frequency, thus providing more relief in the force-limited frequency range. Addi¬ 
tionally, by lowering the fundamental frequencies of the fully-integrated system, isolators 
reduce the transmissibility to the components, whose resonant modes are subsequently in 
the isolation region. It is also known that softer springs, or reducing the stiffness, are more 
effective in the isolation region. The expected increase in damping from the isolators was 
also desired to minimize the peak response at the fundamental frequencies. Therefore, the 
NPSCuL designs were updated to accommodate two different isolators, and vibration tests 
were performed with and without FLVT on the various configurations to demonstrate the 
combined use of FLVT and COTS isolators as a feasible vibration reduction method for 
small satellites. 


3.1 Design Updates to NPSCuL 

The isogrid NPSCuL test results showed that changes to the adapter-ring-to-baseplate in¬ 
terface were the most effective in changing the fundamental frequencies of the system; 
therefore, NPSCuL was re-designed to incoporate isolators at the joint between the adapter 
ring and baseplate. Both this joint and the interface between the P-PODs and the NPSCuL 
wall were considered. However, incorporating isolators at this joint was thought to be the 
most global solution for other ESPA payloads, and existing isolation systems for larger 
spacecraft are implemented at similar locations. 

3.1.1 Design Requirements 

The design changes for incorporating isolators on NPSCuL were driven by the Atlas V 
launch vehicle requirements because it is the most likely launch vehicle that will fly an 
NPSCuL in the future. The pertinent requirements for this design change are as follows [2]: 


41 




• AP static envelope must be within the auxiliary payload (AP) volume 

• AP mechanical interface must be compatible with the ABC plate 

• AP mass properties must be within the interface control document (ICD) range 

• AP must be able to withstand thermal, static loads, and random vibration environ¬ 
ments 

• AP fundamental frequency must be above 35 Hz 

• Mission-specific coupled loads analyses (CLAs) performed by ULA must show no 
negative impacts to the LV and the primary spacecraft 

Other general design considerations included ease of assembly and access to the ABC 
interface fasteners. The following requirements are taken from the ABC User’s Guide and 
do not necessarily reflect the mission-specific requirements for past NPSCuL missions. 

3.1.1.1 Static Envelope 

The non-separating AP envelope for the ABC secondary payload launch system, shown 
in Figure 3.1, had to be maintained for any changes made to NPSCuL. The envelope was 
defined with a tight tolerance (within 2.54 cm, or 1") on the outer dimensions of an NPSCuL 
integrated with eight P-PODs, so there was not much room for expansion in any given 
direction. 
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Figure 3.1: Non-separating AP Envelope. 

Source: Aft Bulkhead Carrier Auxiliary Payload User’s Guide, United 
Launch Alliance Std., May 2014. 


3.1.1.2 Mechanical Interface to Launch Vehicle 

The ABC plate is mounted to the aft end of the Centaur stage with struts (see Figure 1.6). 
The ABC plate interface, shown in Figure 3.2, consists of 24 fasteners on a 38.1 cm (15") 
diameter circle. This is derived from the ESPA standard. 
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+ Y AP 



Figure 3.2: ABC Interface. 

Source: Aft Bulkhead Carrier Auxiliary Payload User’s Guide, United 
Launch Alliance Std., May 2014. 


3.1.1.3 Mass Properties 

The mass properties requirements are provided as a range for both the mass and center 
of gravity (CG). The mass requirement is 65.8 kg ± 11.4 kg (145 ± 25 Ibm), the Z CG 
location must be 24.4 cm ± 1.91 cm (9.6" ± 0.75"), and the X and Y CG offsets must be 
0 cm ± 1.27 cm (0" ± 0.5") [2], However, the mass requirement is often mission-specific 
because it depends on the mass budget and loads analyses results. 

3.1.1.4 Launch Environments 

3.1.1.4.1 Thermal Environment Between the pre-launch and flight thermal environ¬ 
ments, the minimum and maximum expected temperature range is between -46 °C and 
38 °C (-50 and 100 °F). However, the thermal vacuum qualification test range from of -34 
to 71 °C (-29 to 160 °F) were used as a more conservative design requirement to cover any 
extreme temperatures that may occur during flight [2], [12]. 


44 










3.1.1.4.2 Static Loads The launch vehicle design load factors are specified to be 7 Gs 
and 5 Gs in the axial and lateral directions, respectively, and are to be applied at the CG of 
the AP simultaneously. It is noted that unlike most primary spacecraft that are subjected to 
an axial load that is mainly in compression, the axial static load is mainly in tension due to 
the orientation of the AP on the ABC plate. 


3.1.1.4.3 Random Vibration The random vibration environment defined at the ABC in¬ 
terface is shown in Figure 3.3; this is the maximum predicted environment (MPE) by the 
launch vehicle provider for this location. 



Frequency (Hz) 


Figure 3.3: Aft Bulkhead Carrier (ABC) Vibration Requirement. 

Source: Aft Bulkhead Carrier Auxiliary Payload User’s Guide, United 
Launch Alliance Std., May 2014. 


3.1.1.5 Fundamental Frequency 

The launch vehicle also commonly requires that the fundamental frequency of the AP be 
above 35 Hz; however, this general guideline was relaxed for the purposes of choosing an 
isolator, and the impacts of having an isolated AP will require further evaluation by ULA. 
Isolators typically lower the fundamental frequencies of the system, but an effort was made 
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to choose isolators that had higher fundamental frequencies to mitigate the risk of coupling 
between NPSCuL and the Centaur aft bulkhead components. 

3.1.1.6 Mission-Specific CLA 

The CLAs are performed on the entire launch vehicle model, which includes the primary 
spacecraft and AP models, to assess potential impacts to the LV and the primary spacecraft 
for different launch events [2], The CLAs also confirm that the design load factors are 
adequate. A finite element model of the AP is delivered to ULA for these analyses; this has 
been done for all previous missions and will also be done for any future flights that include 
isolators on NPSCuL. 

3.1.2 Isolator Selection 

The two isolators that were chosen for feasibility testing were the LORD conical broad- 
temperature-range (BTR) silicone isolators (P/N AM-009-14) and the Barry cupmount hi- 
damp silicone elastomer isolators (P/N NC-1035-T4). Due to the size and static load ratings 
of the isolators, it was determined that eight isolators were sufficient for both types. 

3.1.2.1 Conical Isolator Description 

The conical isolator shown in Figure 3.4 was initially chosen for feasibility testing due to 
its low profile of 1.80 cm (0.71"), low mass of 81.6 g (0.180 Ibm), and wide operating 
temperature range of -56 to 149 °C (-65 to 300 ° F). It had the highest maximum static load 
per mount of 111 N (25.0 Ibf) and the highest axial natural frequency of 42 Hz compared 
to other similar isolators. Additionally, it was capable of handling the random vibration 
environment shown in Figure 3.3. However, it is generally recommended for use in com¬ 
pression only; its use in tension is to be assessed on a case-by-case basis. This type of 
isolator is also not fail-safe, meaning that if the elastomer fails, there is no other retention 
mechanism to prevent the isolator from falling apart [30]. 
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Figure 3.4: Conical Isolator. 

Source: “LORD aero catalog,” 2011. [Online], Available: 

http://www.lord.com/productsand-solutions/vibration-and-motion- 
control/aerospace-catalog.xml 


3.1.2.2 Cupmount Isolator Description 

The cupmount isolator shown in Figure 3.5 was chosen for its fail-safe and all-attitude 
features. Not only does it have a retention mechanism if the elastomer fails, the isolator can 
support static weight in any direction or orientation; it is not restricted to use in compression 
only [31]. Similar to the conical isolator, the axial natural frequency is 40 Hz, it has a broad 
operating temperature range of -55 to 149 °C (-67 to 300 ° F), and it can support up to 311 
N (70 Ibf) of static load in vibration or 107 to 169 N (24 to 38 Ibf) shock. The cupmount 
isolator weighs 170 g (0.375 Ibm) and is heavier than the conical isolator [32], 



Figure 3.5: Cupmount Isolator. 

Source: “Barry Controls - Core Spin,” June 2015. [Online]. Available: 
http://www. hutchinsonai.com/uploads/tech/Core%20spin. pdf 


Due to its fail-safe and all-attitude features, the cupmount isolator is more suitable for 
NPSCuL on the ABC plate; however, both isolators will be discussed due their potential 
use in other applications that may not have the tension-loading requirement. 
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3.1.3 Design Change to NPSCuL 

The adapter ring and baseplate, shown in Figure 3.6, were modified to accommodate isola¬ 
tors on NPSCuL and to keep the overall dimensions within the static envelope. 



Figure 3.6: Baseline NPSCuL Adapter Ring and Baseplate 


The adapter ring was re-designed to be a flat ring with two different thicknesses. The 
thinner portions of the ring were kept at a 0.432 cm (0.170") thickness, the same thickness 
of the -Z flange of the baseline adapter ring, to maintain the same thread engagement length 
of the interface fasteners between the ring and the ABC plate. The thicker portions of 
the ring were adjusted to accommodate an adequate amount of thread engagement for the 
fasteners holding the isolators to the ring, as well as any nuts or fasteners that needed to be 
installed on the underside of the isolator. The inner and outer diameters of the re-designed 
ring were also adjusted to accommodate the footprint of the isolator. A separate ring was 
designed for each isolator; the three designs are shown in Figure 3.7. The re-designed 
adapter rings, with the isolators installed, are shown in Figure 3.8. In both adapter ring 
designs, only eight of the 24 available ABC interface fasteners are used. 
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(a) Baseline Adapter Ring 


(b) Adapter Ring for Conical 
Isolators 


(c) Adapter Ring for Cupmount 
Isolators 


Figure 3.7: Adapter Ring Designs 



(a) Conical Isolators on Adapter Ring 


(b) Cupmount Isolators on Adapter Ring 


Figure 3.8: Isolators on Re-designed Adapter Rings 


The 24 threaded holes in the baseplate for the ring-to-baseplate fasteners were replaced 
with eight counterbore holes; the size of the counterbore was dependent on the size of 
the isolator-to-baseplate fasteners; the conical isolators required 9.525-mm-diameter (3/8") 
fasteners, and the cupmount isolators required 6.35-mm-diameter (1/4") fasteners. The 
modified baseplate is shown in Figure 3.9. 
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Figure 3.9: Modified Baseplate and Adapter Ring, Cupmount Isolators 


The resulting reductions in overall height are 2.29 cm (0.9") and 1.78 cm (0.7") with the 
conical and cupmount isolators. The adapter ring masses for the conical and cupmount 
isolators are 1.66 and 1.04 kg (3.67 and 2.30 Ibm), respectively. Although the re-designed 
adapter rings are about the same mass or lighter than the baseline adapter ring (weighing 
1.54 kg, or 3.40 lbs), the combined mass with the isolators increased the overall mass of the 
NPSCuL structure by 1.04 and 2.27 kg (2.3 and 2.5 Ibm) with the conical and cupmount 
isolators, respectively. 

3.2 Test Set-Up and Configurations 

3.2.1 FLVT Test Set-Up 

The FLVT set-up shown in Figure 3.10 depicts four force transducers are sandwiched be¬ 
tween the shaker adapter plate and a second plate, or the FLVT adapter plate, that interfaces 
with the test article. As recommended in [23], the mass of the FLVT adapter plate was de¬ 
signed to be no more than 10% of the test article mass, or 8.62 kg (19 Ibm). This set-up 
allows the force transducers, which must be pre-loaded [33], to measure the force at the 
mounting interface. The force-limited set-up was used in all test configurations (with and 
without force-limiting) and for all sine and random test inputs; this was done to maintain a 
consistent test boundary condition. 
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Figure 3.10: FLVT Set-Up, Bare Fixture 


3.2.1.1 Force Limit Parameters 

When using the semi-empirical method, the physical total mass (Mq), the acceleration spec¬ 
tral density (Saa), and the roll-off (fundamental) frequency (/o) are either pre-defined or 
properties of the system. Saa is a given mission requirement, Mq is measured before the 
test, and /o is also measured from the sine sweeps performed before each random test. The 
term and the cut-off frequency, f, are determined by recommendations in [23] or rules 
of thumb based on existing data. For NPSCuL, the recommended value of 2 and cut-off 
frequency of 500 Hz were used in calculating the force limit [23]. It is also recommended 
that the maximum roll-off slope should not exceed 9 dB/Oct, but this recommendation was 
not followed for NPSCuL because the resulting notched ABC acceleration input was de¬ 
termined to be conservative compared to the flight environment, despite a steeper roll-off 
slope. The Matlab code used for the force limit calculation can be found in Appendix E. 

3.2.2 Test Article Description 

For all test configurations, eight P-POD mass models (P2M2s) and the SAD EDU were 
installed on the main NPSCuL structure; this is the typical NPSCuL qualification test con- 
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figuration and is shown in Figure 3.11. Although using P-PODs would have been a more 
flight-like representation of the system, the P2M2s are simple structures than can be eas¬ 
ily modeled and have predictable dynamic behavior; the P-PODs are much more complex 
with variable behavior above 500 Hz due to the internal constraints on the CubeSats and the 
variation in CubeSat dynamic properties. Past data also shows that the acceleration mea¬ 
sured at the P-POD to NPSCuL interface using P2M2s results in a worst-case response; 
this response was measured just below the upper-left fastener for each P-POD location. 



Figure 3.11: Typical Qualification Test Set-Up 


The test configurations are summarized in Table Table 3.1. The conical isolators and mod¬ 
ified adapter ring were installed on both the baseline and isogrid NPSCuL structures (Con¬ 
figurations 2 and 5) for the initial round of feasibility testing. This testing was repeated 
with the cupmount isolators installed on the baseline structure (Configuration 3), and the 
qual test was only performed on this configuration because the cupmount isolator is more 
suitable for NPSCuL on the ABC plate and had the most favorable results observed during 
feasibility testing. The cupmount isolators were not tested with the isogrid structure due to 
behavior observed during the feasibility testing performed with the conical isolators. The 
unmodified baseline and isogrid NPSCuL designs (Configuration 1 and 4) were subjected 
to the same inputs as the other configurations to serve as the comparison point for the 
isolated designs. The baseline NPSCuL was successfully qualification-tested for previous 
missions in June 2011 and September 2012. 
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Table 3.1: Test Configurations 


Configuration 

NPSCuL Structure 

Isolators 

Feasibility Testing 

Qual Test 

1 

Baseline 

N/A 

Y 

Y 

2 

Baseline 

Conical 

Y 

N 

3 

Baseline 

Cupmount 

Y 

Y 

4 

Isogrid 

N/A 

Y 

N 

5 

Isogrid 

Conical 

Y 

N 


3.2.3 Sine and Random Inputs 

Each test configuration was subjected to sine sweeps before and after each random test; this 
cycle of tests was repeated for several random inputs in each axis. The general test flow is 
shown in Figure 3.12. 


Repeated For Each Random Input 



Figure 3.12: General Test Flow for Each Test Configuration 


The sine sweeps were performed at 0.25 or 0.50 G from 20-2000 Hz. In some cases, the 
starting frequency was adjusted down to 5 or 10 Hz to capture lower modes. The pre- and 
post-random sine sweep data were compared at each accelerometer location and the force 
measurement to check for frequency and amplitude shifts of the fundamental modes, which 
are typical indicators of a major change to the system. The sine sweep data was also used 
to experimentally determine the frequency and damping of each fundamental mode, and 
to calculate the apparent mass as a function of frequency. The apparent mass is used as 
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an indication of the transmissibility of the system; in the rigid-body region, it is used as 
a sanity check to ensure that the force transducers have not lost any pre-load [33]. The 
apparent mass is calculated by dividing the measured force by the acceleration control at 
each frequency: 


m((o) 


F{(o) 

a{(D) 


The damping was determined using the half-power bandwidth method [34]: 


(3.1) 



where ^ is the damping ratio, f„ is the natural frequency of interest, and A/ is the frequency 
bandwidth between the -3 dB points, or half-power points, on the FRF. 

The ABC and Operationally Responsive Space (ORS) random inputs, shown in Fig¬ 
ure 3.13, were used in the feasibility testing to investigate the effects of FLVT on the 
isolated and un-isolated spacecraft; the test configurations were subjected to these ran¬ 
dom inputs with and without force limiting. The results from the random vibration tests 
for the ABC profile with force limiting were used for the isolator performance assessment 
because it was a driving requirement for the re-design of NPSCuL. The responses from 
the random vibration tests for the ABC profile without force limiting were used to assess 
the isolators’ effect on NPSCuL without the expected benefit of FLVT. The ORS random 
vibration profile was used to investigate the effects of the isolators and FLVT for a random 
vibration input that has vast differences in frequency content compared to the ABC profile. 
However, it is defined at the LV to P-POD interface and is not realistic for a spacecraft such 
as NPSCuL; the profile was only used due to its high-frequency content. Additionally, the 
target level was adjusted to MPE -6 dB for the ORS profile to avoid over-driving the shaker 
above 400 Hz, where there is a significant amount of energy required to achieve the levels 
for a 90.7 kg (200 Ibm) spacecraft instead of the intended 9.07 kg (20 Ibm) spacecraft. 
The ABC profile was representative of other vibration profiles for similar spacecraft, so no 
other profiles were used; this comparison is shown in Figure 3.14. The summary of inputs. 
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in order of when each configuration was tested, is shown in Table 3.2 and Table 3.3. The 
sine and random vibration tests were performed for each configuration using the test flow 
shown in Figure 3.12. The random input levels used are listed in relative magnitude to the 
maximum predicted environment (MPE) of each vibration profile. 



—ABC MPE (7.6 Grms) — ORS-3, Wafer 2 (P-POD) 15.6 Grms 

Figure 3.13: ABC and ORS Vibration Profiles. 

Adapted from Aft Bulkhead Carrier Auxiliary Payload User’s Guide, United 
Launch Alliance Std., May 2014 and “Operationally Responsive Space 
CubeSat-to-P-POD Interface Control Document, ” Apr. 2013. 
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- ABC MPE (7.6 Grms) - ESPA RUG (MPE), 11.4 Grms 

— NASA GEVS (200 lb payload) 5.6 grms 


Figure 3.14: ABC, ESPA RUG, and GEVS Vibration Profiles. 

Adapted from Aft Bulkhead Carrier Auxiliary Payload User’s Guide, United 
Launch Alliance Std., May 2014; Evolved Expendable Launch Vehicle Sec¬ 
ondary Payload Adapter Rideshare Users Guide (ESPA RUG), DoD Space 
Test Program Std., May 2010; and General Environments Verification Stan¬ 
dard (GEVS), NASA Goddard Space Flight Center Std., Apr. 2013. 


Table 3.2: Sine Test Inputs for All Test Configurations 


Configuration 

NPSCuL 

Structure 

Isolators 

Sine Sweep Frequency Range 

1 

Baseline 

N/A 

0.5 G, 

20-2000 Hz 

4 

Isogrid 

2 

Baseline 

Conical 

5 

Isogrid 

3 

Baseline 

Cupmount 

0.25 G, 

10-2000 Hz (Z) 

5-2000 Hz (X, Y) 
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Table 3.3: Random Test Inputs for All Test Configurations 


Configuration 

NPSCuL 

Structure 

Isolators 

ABC Levels 

ORS Levels 

Unnotched 

FLVT 

Unnotched 

FLVT 

1 

Baseline 

N/A 

MPE 

+0 dB, 

20 sec 

MPE 

+0dB, 

20 sec 

MPE 

-6 dB, 

20 sec 

N/A 

4 

Isogrid 

2 

Baseline 

Conical 

5 

Isogrid 

3 

Baseline 

Cupmount 

MPE 

+6dB, 

3 min 

MPE 

-6dB, 

20 sec 
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CHAPTER 4: 
Experimental Results 


The results from the baseline and isogrid NPSCuL test configurations (Configurations 1 
and 4) were used as comparison point for evaluating the effectiveness of the isolators and 
force limiting. These two configurations had been previously subjected to the ABC input 
with force limiting [19], but the additional testing performed using the remaining three 
random vibration profiles allow for consistent comparisons between the isolated and un¬ 
isolated configurations. Although the cupmount isolator configuration was subjected to the 
higher, qualification-level random vibration test, the responses measured during MPE -i-O 
dB for all configurations were used for comparison and performance evaluation purposes. 

The evaluation metrics were the apparent mass measured at the base of NPSCuL, the en¬ 
veloped acceleration response at the P-POD interface, and the corresponding Grms- The 
apparent mass was used to determine the measured dynamic properties as well as an in¬ 
dication of the response at the P-POD interface. The FRF is inversely proportional to the 
impedance (Eqn. 2.6), so a lower apparent mass indicates that the impedance is high, thus 
resulting in a lower response at the P-POD interface. The enveloped acceleration data are 
the max ASD values of every instrumented P-POD location in the specified axis; it includes 
the cross-axis response as well. For example, the X-axis envelope would be the maximum 
ASD values measured in the X-axis, at the interface of all eight P-PODs, during the X, Y, 
and Z-axis tests. Enveloping the response by axis is the most straight-forward method of 
assessing the isolators’ performance and is independent of the indivudual P-PODs’ orien¬ 
tations on NPSCuL. The data for the individual P-POD responses (in Grms) can be found 
in Appendix C. 

4.1 Conical Isolator Results 

4.1.1 Conical Isolators on Baseline NPSCuL 

The apparent mass plots of the baseline NPSCuL with and without conical isolators (Con¬ 
figurations 1 and 2) are shown in Figure 4.1, and the fundamental frequencies measured 
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from these sine sweeps are shown in Table 4.1. With the isolators (Configuration 2), there 
are two distinct modes in each lateral axis (X and Y) that are much lower in amplitude com¬ 
pared to the fundamental mode in the baseline NPSCuL. These modes correspond to the 
rocking and shear modes in the lateral directions. This is typical behavior for base-mounted 
isolators on a box-like structure; many isolated avionics boxes on launch vehicles exhibit 
similar dynamic behavior [35]. In the axial direction (Z), there are decreases in both the 
frequency and the amplitude of the fundamental mode. The decrease in frequency, which 
was expected, is due to the stiffness properties of the isolators. The reduction in amplitude 
of all of the fundamental modes can be attributed to the increased damping provided by 
the isolators. As expected, the isolation region in each axis starts at approximately V^cOn, 
where (On is the second fundamental frequency in the X and Y axes and the first funda¬ 
mental frequency in the Z axis. In general, the apparent mass above 500 Hz is very low, 
indicating that the response at the P-POD interface will be lower as well. The peak at 1560 
Hz in the Z-axis corresponds to an FLVT plate mode that is excited in the axial direction. 
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Frequency, Hz 



Figure 4.1: Apparent Mass, Conical Isolators on Baseline NPSCuL 
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Table 4.1: Measured Frequencies from Apparent Mass, Conical Isolators 
on NPSCuL 


Mode 

Measured Frequencies (Hz) 

Configuration 

Configuration 


1 

2 

1 

(Rocking, Y) 

49.7 

19.8 

2 

(Rocking, X) 

48.5 

20.4 

3 

(Axial, Z) 

117.0 

56.3 


The enveloped acceleration responses at the P-POD interface in each test axis for both 
configurations are shown in Figure 4.2; these responses were measured during the random 
vibration test for the ABC profile with force limiting. The amplitudes of the responses 
up to 200 Hz are comparable in all axes, and the peaks are shifted due to the differences 
in fundamental frequencies. However, the shift in frequency above 200 Hz is negligible, 
and the response consistently drops by two orders of magnitude in all axes. This indicates 
that COTS isolators are effective in reducing the P-POD response response above 500 Hz, 
which provides the most benefit to the CubeSats. 
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Config 1 —Config 2 200 Hz 


Frequency (Hz) 

Config 1 —Config 2 200 Hz 




(a) X-Axis Envelope (b) Y-Axis Envelope 



Config 1 — Config 2 200 Hz 


(c) Z-Axis Envelope 

Figure 4.2: Acceleration Responses at P-POD Interface Due to ABC Input 
with FLVT, Baseline NPSCuL with Conical Isolators 


The Grms from 20 to 2000 Hz and from 500 to 2000 Hz for both configurations, shown 
in Table 4.2 and Table 4.3, confirm that the high-frequency content (above 500 Hz) is 
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reduced by up to 98% when isolators are used. These values correspond to the data shown 
in Figure 4.2. As shown in Figure 4.3, the high-frequency content accounts for up to 50% 
of the Grms on the baseline NPSCuL; however, the high-frequency content is no more than 
10% of the total Grms when conical isolators are incorporated. 

Table 4.2: Grms Reduction at P-POD Interface, 20-2000 Hz — Conical 
Isolators 


Test 

Axis 

Grms 

Percent 

Difference, % 

Configuration 1 

Configuration 2 

X 

15.9 

6.15 

61.2 

Y 

16.7 

5.47 

67.3 

Z 

8.41 

4.36 

48.1 


Table 4.3: Grms Reduction at P-POD Interface, 500-2000 Hz — Conical 
Isolators 


Test 

Axis 

Grms 

Percent 

Difference, % 

Configuration 1 

Configuration 2 

X 

7.73 

0.601 

92.2 

Y 

7.46 

0.487 

93.5 

Z 

2.33 

0.169 

92.7 
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Figure 4.3: Reduction in High Frequency Content at P-POD Interface, 
Conical Isolators 


4.1.1.1 Effect of Fully Constrained Conical Isolators on NPSCuL 

Initially, the conical isolators were mistakenly installed as shown in Figure 4.4; these re¬ 
sults were presented in [1]. In this configuration, the elastomer was fully constrained (in 
compression, tension, and shear) because the fastener that joins the supporting unit (base¬ 
plate) to the isolator was also threaded into the supporting member (adapter ring). This 
prevented the elastomer from deflecting as intended by the manufacturer, resulting in little 
reduction in joint stiffness at each isolator location. 

□ Elastomer 

E3 Metal Core 

Figure 4.4: Apparent Mass, Conical Isolators on Baseline NPSCuL 



Figure 4.5 shows the apparent mass plots of this configuration with fully-constrained con¬ 
ical isolators in comparison with those of the baseline NPSCuL. Although the first funda- 
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mental frequency in the X and Y axes dropped by approximately 15 Hz, there is no signi¬ 
fication reduction in overall amplitude. The reduction in overall joint stiffness at the ring- 
to-baseplate interface, not the properties of the isolators, are attributed to these changes. 
Because only eight of the original 24 fasteners of the LV interface are used, and the contact 
surface area between the baseplate and the adapter ring is significantly less, the overall joint 
stiffness at this interface is lower than that of the baseline NPSCuL. In the Z-axis, both the 
apparent mass amplitude and the fundamental frequency are lower, which indicates that 
the reduction in overall joint stiffness may also be a result of the isolator properties in this 
direction. The apparent mass plots indicate that the conical isolators are not effective in 
the lateral directions when fully constrained. However, this constraint method is a possi¬ 
ble alternative in the event that an all-attitude isolator cannot be used in a static tension 
application and the significant reduction in isolation performance is acceptable. 
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Figure 4.5: Apparent Mass, Fully Constrained Conical Isolators on Baseline 
NPSCuL 


4.1.2 Conical Isolators on Isogrid NPSCuL 

The apparent mass plots shown in Figure 4.6 exhibit very similar responses between the 
baseline and isogrid NPSCuL when isolators are ineorporated. The mass and CG of 
NPSCuL are designed to be the same with both designs, so it is likely that there are similar 
responses because the unibase on the isgorid NPSCuL was replaced with a separate ring 
and baseplate that is in the same configuration as the baseline NPSCuL. 
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Figure 4.6: Apparent Mass, Conical Isolators on Baseline and Isogrid 
NPSCuL 


Due to the comparable performance between the two configurations and the additional 
work that would be required to use the isgorid structure for a future flight (i.e. accom¬ 
modating the structure for ground support equipment (GSE), additional load and environ¬ 
mental testing), no further testing was performed on the isogrid NPSCuL, and it was not 
considered for future missions. 
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4.2 Cupmount Isolator Results 

The apparent mass comparison plots for the baseline NPSCuL with and without the cup- 
mount isolators (Configurations 1 and 3) are shown in Figure 4.7, and the corresponding 
fundamental frequencies measured from the sine sweeps are shown in Table 4.4, along with 
the frequencies measured for the baseline NPSCuL with conical isolators. Similar to the 
results with the conical isolators, there are two distinct modes in each lateral axis (X and 
Y) that are much lower in amplitude compared to the fundamental mode in the baseline 
NPSCuL, and the apparent mass in the isolation region (above is greatly reduced in 

all axes. 



10 "^ 


10 ^ 10 ^ 10 ® 

Frequency, Hz 




Figure 4.7: Apparent Mass, Cupmount Isolators on Baseline NPSCuL 
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The fundamental frequencies in the lateral directions are also below 20 Hz, where the 
ABC random input starts. The ASD values are not defined below 20 Hz, so the roll-off 
frequency for the testing performed on Configuration 3 with force-limiting is 20 Hz instead 
of the fundamental frequencies shown in Table 4.4. 

Table 4.4: Measured Frequencies from Apparent Mass, Cupmount Isolators 
on NPSCuL 


Mode 

Measured Frequencies (Hz) 

Configuration 

Configuration 

Configuration 


1 

2 

3 

1 

(Rocking, Y) 

49.7 

19.8 

12.4 

2 

(Rocking, X) 

48.5 

20.4 

12.5 

3 

(Axial, Z) 

117.0 

56.3 

39.5 


Compared to the results of the conical isolators, the cupmount isolators are more effective 
due to the lower fundamental frequencies, as shown in Figure 4.8. Although the perfor¬ 
mance is similar above 500 Hz in the X- and Y-axes and 300 Hz in the Z-axis, the isolation 
region of the cupmount isolators starts at a lower frequency in the lateral directions, thus 
indicating a reduction in response over a broader frequency range. 


70 





Frequency, Hz 



-— Config 2 - Config 3 - Start of Isolation Region, Config 2-Start of Isolation Region, Config 3 


Figure 4.8: Apparent Mass, Conical and Cupmount Isolators on Baseline 
NPSCuL 


The enveloped acceleration responses at the P-POD interface in each test axis for the base¬ 
line NPSCuL with and without cupmount isolators are shown in Figure 4.9; these responses 
were measured during the random vibration test for the ABC profile with force limiting. 
The responses at the peak frequencies up to 300 Hz are comparable to the baseline NPSCuL 
in all axes, and the response consistently drops by two orders of magnitude (20 dB) in all 
axes above 300 Hz in the X- and Y-axes and 70 Hz in the Z-axis. It is noted that the re¬ 
sponses are lower below 300 Hz as well, except at the peaks. This is most likely due to the 
significant apparent mass of the second mode, so there is some reduction in amplitude in 
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the lower frequencies, but full isolation does not occur until after this second mode. Similar 
to the conical isolators, the cupmount isolators also effectively reduce the response above 
500 Hz. 




(a) X-Axis Envelope (b) Y-Axis Envelope 



(c) Z-Axis Envelope 

Figure 4.9: Acceleration Responses at P-POD Interface Due to ABC Input 
with FLVT, Baseline NPSCuL with Cupmount Isolators 
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The P-POD interface responses with the cupmount isolators (Configuration 3) were also 
compared to those with the conical isolators (Configuration 2); these plots are shown in 
Figure 4.10. Although the apparent mass above 300 Hz is the same with both isolators, the 
cupmount isolators reduce the response at the P-POD interface even further over the entire 
frequency range from 20 to 2000 Hz. The resulting Grms is also lower by an additional 
32% overall and 57% above 500 Hz, as shown in Table 4.5 and Table 4.6. Figure 4.11 
highlights the reduction in high-frequency content. These results confirm that when the 
isolation region starts at a lower frequency, the isolator is more effective. This ultimately 
results in a better vibration environment at the P-POD interface. 
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ASD (g^/Hz) 




(a) X-Axis Envelope (b) Y-Axis Envelope 



(c) Z-Axis Envelope 

Figure 4.10: Acceleration Responses at P-POD Interface Due to ABC Input 
with FLVT, Baseline NPSCuL with Conical and Cupmount Isolators 
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Table 4.5: Grms Reduction at P-POD Interface, 20-2000 Hz — Cupmount 
Isolators 


Test 

Axis 

Grms 

Percent 

Difference, 

% 


Configuration 

Configuration 

Configuration 

Config. 

Config. 


1 

2 

3 

1 vs 3 

2 vs 3 

X 

15.9 

6.15 

3.76 

76.3 

38.8 

Y 

16.7 

5.47 

3.72 

77.7 

32.0 

Z 

8.41 

4.36 

2.96 

64.9 

32.3 


Table 4.6: Grms Reduction at P-POD Interface, 500-2000 Hz - Cupmount 
Isolators 


Test 

Axis 

Grms 

Percent 

Difference, 

% 


Configuration 

Configuration 

Configuration 

Config. 

Config. 


1 

2 

3 

1 vs 3 

2 vs 3 

X 

7.73 

0.601 

0.257 

96.7 

57.3 

Y 

7.46 

0.487 

0.240 

96.8 

50.6 

Z 

2.33 

0.169 

0.077 

96.7 

54.7 
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Test Axis 


□Config 1 BConfig 2BConfig 3 


Figure 4.11: Reduction in High Frequency Content at P-POD Interface, 
Cupmount Isolators 


4.2.1 Isolator Characterization - Test Like You Fly 

Due to the orientation of NPSCuL on the ABC plate, the driving static loading condition is 
in tension; this is not a typical loading condition for spacecraft. The practice of testing in 
a flight-like condition, or "test like you fly," is a primary concern in the aerospace industry 
because it is often difficult to replicate flight conditions on the ground for various envi¬ 
ronments. The vibration tests performed on the ground are usually done in a compressive 
static loading condition, but it was not practical to replicate the tension loading condition 
for NPSCuL. Although the cupmount isolators are all-attitude, a separate characterization 
test was performed to verify that the isolator performance was comparable in both com¬ 
pression and tension, thus alleviating the concern that the typical vibration test set-up may 
not be valid for NPSCuL on the ABC. 

The test set-up for the tension and compression cases in the axial (Z) direction are shown 
in Figure 4.12 and Figure 4.13, in which a P2M2 is supported by four cupmount isolators. 
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The tension and compression configurations were tested in the axial (Z-axis) and worst- 
case lateral (Y-axis) directions. The isolators and mass model were mounted at the same 
height for both cases to maintain a consistent boundary condition at the isolator input loca¬ 
tion. The input accelerometer for the Z-axis test was located in the middle of the isolator 
mounting plate, and the input accelerometer for the lateral test was located on the shaker 
adapter plate for better control. The test set-ups for the lateral direction are shown in Fig¬ 
ure 4.14 and Figure 4.15. In both set-ups, the accelerometers were located on both the 
isolator mounting plate and the mass model interface plate to measure the response across 
the isolators. Following the test flow in Figure 3.12, sine sweeps and a 3-minute random 
vibration test at MPE -i-O dB of the ABC input levels were performed for all set-ups shown. 



,^^Shaker 
Adapter Plate 


• Mass Model 
Interface Plate 


^P-FOD 
Mass Model 
(P2M2) 


Isolator 

Mounting 

Plate 


Figure 4.12: Z-Axis Test, Tension 
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Figure 4.13: Z-Axis Test, Compression 



Figure 4.14: Lateral Test, Tension 
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Figure 4.15: Lateral Test, Compression 


As shown in Figure 4.16 and Figure 4.17, the acceleration responses measured on the mass 
model interface plate in both compression and tension tests are nearly identical for both 
axial and lateral test configurations. This indicates that the test configuration for NPSCuL 
represents a test-like-you-fly configuration with respect to dynamic loads, despite the dif¬ 
ference in static the static loading conditions between the flight and test configurations. 
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[g'/Hz] 



Figure 4.16: Lateral Random Vibe Test, Compression and Tension at MPE 
+0 dB 



20 100 1000 2000 


[Hz] 


Figure 4.17: Lateral Random Vibe Test, Compression and Tension at MPE 
+0 dB 
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4.2.2 Break-In Testing 

As recommended by the vendor, the cupmount isolators were subjected to a "break-in" test 
to relieve the molded-in stresses in the elastomer. Unlike the conical isolators, the elastomer 
is not bonded to the metal inner core or outer cup; therefore, there is some residual stress 
from molding the elastomer in the cupmount isolator. The break-in test can be performed 
by deflecting the isolator in the axial direction to at least three times the published deflection 
limits; this can be done by using a hand-operated arbor press instrumented with a load cell 
and deflection gage. Alternatively, the break-in test can also be performed by running an 
extra operational-level vibration test prior to testing [36]. The break-in test performed for 
the cupmount isolators on NPSCuL was a 20-second duration qualification-level (MPE -1-6 
dB) random vibration test in the Z-axis. The pre- and post-random sine sweeps, plotted in 
Figure 4.18, show a shift in fundamental frequency of 20 Hz. The subsequent sine sweeps 
did not exhibit any additional shifts, thus indicating that the isolators were properly broken 
in. 



- PreBreak-ln,Down - PostBreak-ln(PreABC),Down 


Figure 4.18: Pre- and Post-random Sine Sweep Comparison, Z-Axis Break- 
In Test 
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4.2.3 Qualification Testing 

Based on the favorable results of this configuration during feasibility testing and the cup- 
mount isolator’s all-attitude, fail-safe features, the baseline NPSCuL with cupmount iso¬ 
lators was subjected to a qualification-level test in which the random vibration test was 
performed to MPE +6 dB of the ABC profile for a duration of 3 minutes. Performing 
the qualification test provided confidence that incorporating these isolators on NPSCuL 
would be a viable option for future flights. The fasteners that go through the baseplate 
and into the threaded core of the isolators (see Figure 3.5) must be torqued appropriately 
for the expected dynamic loading condition. It was determined that a torque of 125 in-lbs 
was required, but this was not achievable without a means of restraining the threaded core; 
anything higher than approximately 60 in-lbs caused the threaded core to spin. The vendor- 
recommended torque value of 50 in-lbs was based on a lower grade bolt (Grade 2 instead 
of A286 stainless steel) and was designed to prevent the threaded core from spinning. Due 
to a lack of experience with these isolators, the vendor-recommended torque value was 
applied without any secondary back-out prevention in the initial qualification test, and six 
of the eight isolators were damaged due to insufficient torque and subsequent gapping; the 
damage is shown in Figure 4.19. 



Visible 

damage 

to 

threaded 

core 



Figure 4.19: Damaged Isolators Due to InsufFicient Torque 


No foreign-object-debris (FOD) resulted from the damage. The post-test disassembly was 
performed after the structure had been qualification tested in all three axes; at this point, 
it was discovered that the fasteners had backed out of, but were not able to full disengage 
from, the core due to the low clearance between the P2M2s and the baseplate. Two of the 
isolators were still intact, and there was enough thread engagement to keep the structure 
together. Additionally, there was no clear indication in the test data that indicated a failure 
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had occurred. The initial drop in fundamental frequency was attributed only to the break- 
in process, not any loss in pre-load, and the subsequent shifts in fundamental frequency 
were within the typical 10% range. It is noted that time histories were not saved due to the 
limitations of the data acquisition system; all data captured was in the frequency domain. 
The availability of the time histories might have aided in determining when the damage 
oceurred. 

The qualification test was successfully repeated with a new set of isolators; the fasteners 
were torqued to 14.1 N-m (125 in-lbf) and thread-lock was applied to each of the eight fas¬ 
teners for secondary backout prevention. To aehieve the proper torque value, the threaded 
eore was restrained with another fastener on the underside of the isolator; this allowed 
equal and opposite moments to be applied on both sides. 

4.2.4 Temperature Effects 

The isolator temperatures were measured to alleviate the eoneern that the frequencies would 
shift significantly due to a large change in temperature, induce coupling, and cause failures 
to the launch vehicle or primary payload. Thermocouples were attached to six of the eight 
isolators during the qualifieation test; the isolator designations are shown in Figure 4.20, 
and the thermocouple taped to Isolator 1 is shown in Figure 4.21. 
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Figure 4.20: Isolator Designation 



Figure 4.21: Thermocouple on Isolator 1 


Temperatures were recorded for all three axes at five-second intervals, as shown in Fig¬ 
ure 4.22 through Figure 4.24. It is noted that no temperature data was captured during 
the Y-axis sine-sweep. The maximum temperature difference observed was 12 °C (22 °F) 
during the Y-axis qual test, from the beginning of MPE -18 dB through the end of MPE 
-1-0 dB of target level. However, the maximum temperature difference observed through 
the 20-second-duration acceptance level (ABC MPE -i-O dB) was only 0.56 °C (1 °F). The 
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temperature difference may have been larger for a full 60-second-duration acceptance test, 
but the data indicates that it likely would not be more than 2.2 °C (5 °F). Also, there was no 
appreciable difference in peak frequencies and amplitudes during the 3-minute qualifica¬ 
tion test despite the temperature increase. The measured fundamental frequencies in each 
axis are shown in Table 4.7; all shifts are within 10%. 

Beginning of Acceptance 
End of Acceptance/Beginning of Proto-Qual 
Beginning of Qual/End of Proto-Qual 
End of Qual 
-Isolator 1 
-Isolator 2 
-Isolator 3 
-Isolator 4 
-Isolator 5 
-Isolator 6 

0 200 400 

Time Elapsed (s) 

Figure 4.22: Isolator Temperatures — MPE 4-6 dB, Z 
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Temperature (°C) Temperature (°C) 



Beginning of Acceptance 
End of Acceptance/Beginning of Proto-Qual 
Beginning of Qual/End of Proto-Qual 
End of Qual 
Start of Sine Sweep 
Isolator 1 
-Isolator 2 
Isolator 3 
-Isolator 4 
-Isolator 5 
-Isolator 6 


Time Elapsed (s) 


Figure 4.23: Isolator Temperatures — MPE +6 dB, X 
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Figure 4.24: Isolator Temperatures — MPE +6 dB, Y 
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Table 4.7: Measured Fundamental Frequencies, Pre- and Post-random Sine 
Sweeps 


Axis 

Mode 1 

Mode 2 

Mode 3 

X 

Pre-Qual 
Frequency (Hz) 

12.3 

63.5 

79.8 

Post-qual 
Frequency (Hz) 

12.5 

65.4 

82.2 

Percent 

Difference (%) 

1.63 

2.99 

3.01 

Y 

Pre-Qual 
Frequency (Hz) 

11.4 

60.3 

84.6 

Post-qual 
Frequency (Hz) 

12.4 

62.8 

85.6 

Percent 

Difference (%) 

8.77 

4.15 

1.18 

Z 

Pre-Qual 
Frequency (Hz) 

43.3 

N/A 

Post-qual 
Frequency (Hz) 

43.3 

Percent 

Difference (%) 

0.00 


The effects of the ambient temperature are also a common concern for elastomeric isolators; 
depending on the material used, the dynamic stiffness can increase by up to nine times at 
cold temperatures. Cold temperatures can also cause the elastomer to be brittle and fail 
under load. However, the cupmount isolators selected for NPSCuL use silicone, which does 
not vary significantly with ambient temperature, especially for the typical ABC thermal 
environment. The stiffening ratio is shown in Figure 4.25. 
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Figure 4.25: Dynamic Modulus Variation of Silicone Due to Temperature. 
Source: “Elastomer dynamic modulus stiffening vs 

temperature,” Oct. 2015. [Online], Available: 

http://www. hutchinsonai.com/uploads/tech/Elastomer_Temp_ Range, pdf 


4.3 Combined Effects of FLVT and Isolation 

The isolators are effective in reducing the response at the P-POD interface on NPSCuL, 
but comparing the results between the tests performed with and without FLVT show that 
the test conservatism is further reduced when FLVT is employed. The envelopes of the 
acceleration responses at the P-POD interface in each test axis for Configurations 1 through 
3 with and without FLVT are shown in Appendix B; these responses were measured during 
the random vibration test for the ABC profile without force limiting. As expected, the high- 
frequency response for the isolated configurations are comparable, similar to the force- 
limited data. The force-limited set-up was not intended or expected to affect the response 
in this frequency range. The corresponding Grms values calculated over the entire test 
frequency range for the unnotched tests are shown in Table 4.8. 
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Table 4.8: X-Axis Grms Envelope of P-POD Interface at MPE +0 dB, 
ABC Profile, Unnotched 


Test 

Axis 

Grms 

Configuration 

1 

Configuration 

2 

Configuration 

3 

X 

29.6 

9.66 

3.43 

Y 

29.9 

10.9 

4.62 

Z 

25.1 

7.35 

3.36 


The Grms values from Table 4.5 and Table 4.8 are plotted in Figure 4.26 to show the 
significant reduction in response due to the combined use of isolators and FLVT. If the fun¬ 
damental modes of the system lie within the random vibration acceleration input frequency 
range of 20-2000 Hz, force limiting is still beneficial in reducing the response at the P- 
POD interface; this is exhibited in the significant reductions for Configurations 1 and 2 
when FLVT is employed. However, FLVT is not as beneficial as the isolators become more 
effective due the shift in fundamental frequency outside of the test frequency range. This 
is most evident in the X and Y axes for Configuration 3, where the fundamental frequency 
is 12 Hz instead of 20 or 50 Hz. 



X Y Z 

Test Axis 

■ Configl BConfigZ ■ ConfigS 

(a) ABC Input, Unnotched 






t 


Test Axis 

I Configl ■ Config2 ■ ConfigS 

(b) ABC Input, FLVT 


Figure 4.26: Overall Grms Envelope at P-POD Interface, MPE -fO dB 
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4.4 Effects of Varying Random Vibration Inputs 

The baseline NPSCuL with and without cupmount isolators were also subjected to the ORS 
random vibration profile without force limiting; only the baseline NPSCuL with cupmount 
isolators was subjected to the force-limited ORS random vibration profile. The Grms values 
calculated over 20 to 2000 Hz corresponding to the tests without force limiting are shown 
in Table 4.9; the Grms values calculated over 20 to 2000 Hz corresponding to the test with 
force limiting are shown in Table 4.10. The envelopes of the acceleration responses at the 
P-POD interface in each test axis due to the ORS random input without FLVT are shown 
in Figure 4.27. 


Table 4.9: Grms Reduction at P-POD Interface, 20-2000 Hz — MPE -6 
dB of ORS Input without Force Limiting 


Test 

Configuration 1, 

Configuration 3, 

Percent Difference 

Axis 

Grms 

Grms 

(%) 

X 

30.3 

3.90 

87.2 

Y 

27.3 

3.68 

86.5 

Z 

12.1 

1.66 

86.2 


Table 4.10: Grms Reduction at P-POD Interface on Configuration 3, 20- 
2000 Hz — MPE -6 dB of ORS Input 


Test 

Axis 

Configuration 3, 

Unnotched 

Grms 

Configuration 3, 

FLVT 

Grms 

Percent Difference 

(%) 

X 

3.90 

3.88 

0.513 

Y 

3.68 

3.73 

1.36 

Z 

1.66 

1.69 

1.81 
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(a) X-Axis Envelope (b) Y-Axis Envelope 



(c) Z-Axis Envelope 

Figure 4.27: Acceleration Response at P-POD Interface Due to ORS Input 
without FLVT, Baseline NPSCuL with and without Cupmount Isolators 


Due to the large amount of energy above 400 Hz in the ORS profile compared to the ABC 
profile, the isolators are very effective in reducing the response at the P-POD interface. 
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However, there is very little energy in the ORS profile below 400 Hz, especially in the 
frequency range of the fundamental modes of the system. Therefore, force limiting is not 
beneficial for this type of input. The Grms values from Table 4.9 and Table 4.10 are plotted 
in Figure 4.28 and Figure 4.29 to illustrate this finding - there is a large reduction in Grms 
between the isolated and baseline configurations, but there is no significant change with 
and without FLVT. 
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25.00 
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I 20.00 

15.00 
10.00 
5.00 
0.00 

X Y Z 

Test Axis 

■ Configl ■ Config 3 



Figure 4.28: Overall Grms Envelope at P-POD Interface, MPE -6 dB of 
ORS Vibration Profile, Unnotched 



Figure 4.29: Overall Grms Envelope at P-POD Interface, MPE -6 dB of 
ORS Vibration Profile, FLVT 
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4.5 Force Limit Example of Isolated NPSCuL 

The plots shown in Figure 2.6 and Figure 2.7 are repeated for the isolated NPSCuL to show 
that the 2DOF assumption in justifying the notched acceleration input is still valid for an 
isolated spacecraft. In the model, the fundamental frequency of the isolated NPSCuL was 
set to 50 Hz, which corresponds to the measured fundamental frequency during the Z-axis 
sine sweep for the baseline NPSCuL with cupmount isolators. The difference between the 
baseline and isolated NPSCuL models is the fundamental frequency, both of which are 
below 240 Hz. In this frequency range, the ASD values of the unnotched acceleration input 
are still relatively high. To properly compare the acceleration response of the ABC plate 
with the force limit used in test, the damping of the isolated NPSCuL model was set to 0.22, 
which was measured from the sine sweep data. The large increase in damping from 0.02 
to 0.22 causes the anti-resonance of the ABC plate in Figure 4.30 to be muted. However, 
Figure 4.31 shows that the 2DOF model is still adequate for predicting the magnitude of 
the notched input near the fundamental frequency of the fixed-base spacecraft. 



—Fixed Base (SDOF), Isolated NPSCuL — 2DOF ABC Plate ,m1 (Isolated NPSCuL)| 


Figure 4.30: Transmissibility of Isolated NPSCuL (SDOF Model) and ABC 
Plate (2DOF Model) 
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Figure 4.31: Predicted and Actual Force-Limited Acceleration Input, Iso¬ 
lated NPSCuL 
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CHAPTER 5: 

Finite Element Model and Analyses 


The finite element models (FEMs) were constructed and used for analysis in NX and NAS- 
TRAN. The baseline NPSCuL FEM that has been delivered for all launch vehicle mission- 
specific CLAs was modified for each isolator configuration. All of the FEMs were corre¬ 
lated with the modal properties measured during testing and provided a sanity check for 
these modal properties. However, only the FEM for the baseline NPSCuL with cupmount 
isolators was used for loads analyses and was delivered to ULA for CLA purposes because 
it is the most likely configuration to fly on a future ABC mission. All analyses were per¬ 
formed with the model fixed at rigid-body element used to tie the nodes representing the 
bolted interface to the ABC plate together. 

5.1 Baseline NPSCuL FEM 

The baseline NPSCuL FEM is shown in Figure 5.1. The primary structure consists of thin- 
shell elements, and the P-POD/P2M2, SAD, and harness mass properties are represented as 
concentrated mass elements that are rigidly connected to the primary structure. All fasten¬ 
ers are modeled as beam elements to facilitate adjusting the joint stiffness for correlation 
purposes and to extract fastener loads from the static load analyses. The modal proper¬ 
ties for this model were correlated to the measured fundamental frequencies from the sine 
sweeps. 
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Figure 5.1: Baseline NPSCuL FEM 


The material properties, joint stiffnesses, and element physical properties used in this model 
are shown in Table 5.1 and Table 5.2. The entire primary structure is made of aluminum 
7075, and all fasteners are A286 stainless steel. The masses associated with the concen¬ 
trated mass elements are listed in Table 5.3, and the overall FEM mass properties are shown 
in Table 5.4. 


Table 5.1: Material Properties 


Material 

Density, 

kg/rn^ 

\lbm/ir?] 

Poisson’s 

Ratio 

Modulus of 

Elasticity (E), GPa 
[psi] 

Aluminum 7075 

2.85E-t03 

[0.103] 

0.33 

68.9 

[1.04E+07] 

A286 Stainless Steel 

7.92E-I-02 

[0.286] 

0.33 

27.6 

[4.00E+06] 
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Table 5.2: Element Physical Properties 


Element 

Property 

Type 

Element Description 

Thickness, 
cm. [in.] 

Non-Structural Mass (NSM), 

kg/n?' [Ibm/in^] 


Angle Brackets 

0.476 

[0.1875] 

0 

Shell 

Side Wall 

0.635 

[0.250] 

2.55 

[3.63E-03] 


Baseplate 

1.27 

[0.500] 

0 


Ring (Upper) 

0.343 

[0.135] 

0 


Ring (Wall) 

0.635 

[0.250] 

0 


Ring (Bottom) 

0.432 

[0.170] 

0 

Bar 

Baseplate-to-Wall Fasteners 

0.635 (dia.) 
[0.250] 

20.7 

[2.94E-02] 


Table 5.3: Concentrated Mass Element Properties 


Element Description 

Mass, kg [Ibm] 

P-POD 1 

8.12 [17.9] 

P-POD 2 

8.12 [17.9] 

P-POD 3 

8.07 [17.8] 

P-POD 4 

8.12 [17.9] 

P-POD 5 

8.12 [17.9] 

P-POD 6 

8.12 [17.9] 

P-POD 7 

8.12 [17.9] 

P-POD 8 

7.94 [17.5] 

SAD 

2.04 [4.50] 

Harnessing 

2.89 [6.37] 
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Table 5.4: FEM Mass Properties Summary 


Property 

Value 

Mass, kg [Ibm] 

85.7 [189] 

Xabc, cm [in] 

-9.96E-03 [-3.92E-03] 

Yabc, cm [in] 

6.68E-03 [2.63E-01] 

Zabc, cm [in] 

26.2 [10.3] 


5.1.1 Fundamental Modes — Frequencies and Mode Shapes 

The fundamental frequencies of the baseline NPSCuL FEM are shown in Table 5.5 along 
with the measured frequencies from the sine sweep test. There is good correlation for the 
rocking modes, but the axial mode of the FEM is not as stiff as the measured frequency. 
This discrepancy was acceptable for the CLAs because the lower fundamental frequency 
in the FEM produces more conservative results when used in the static load analyses. Ad¬ 
ditionally, the frequency range of concern in the CLAs is often cut off around 50 Hz, so the 
large discrepancy was not considered a risk for the LV. The corresponding mode shapes are 
shown in Figure 5.2. 


Table 5.5: Baseline NPSCuL — Fundamental Frequencies 


Mode 

Frequencies (Hz) 

Percent 

Difference (%) 

FEM 

Test 

1 (Rocking, X) 

50.5 

48.5 

4.12 

2 (Rocking, Y) 

50.9 

49.7 

2.41 

3 (Axial, Z) 

90.3 

117.0 

22.8 
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(a) Baseline NPSCuL — Rocking Mode, X-Axis (b) Baseline NPSCuL — Rocking Mode, Y-Axis 



(c) Baseline NPSCuL — Axial "Pogo" Mode, Z- 
Axis 

Figure 5.2: Baseline NPSCuL Fundamental Modes — Mode Shapes 


5.2 Baseline NPSCuL FEM with Conical Isolators 

The baseline NPSCuL FEM was modified to reflect the test configuration with conical 
isolators. As shown in Figure 5.3, the baseline adapter ring was replaced with a flat ring and 
the isolators, which are modeled using thin-shell and spring elements (CBUSH elements in 
Nastran), respectively. The thickness of the raised portions of the ring that accommodate 
the isolators is 1.17 cm (0.460"), and the thickness of the thinner portions of the ring that 


99 








interface with the ABC plate is 0.432 cm (0.170"); these are accounted for in the shell 
properties of the model. The spring stiffness value in all degrees of freedom was 525 kN/m 
(3000 Ibf/in) to correlate the model with the test frequencies. The remaining NPSCuL 
structure and components were not modified. 



Figure 5.3: Flat Ring and Conical Isolator FEM 


5.2.1 Fundamental Modes — Frequencies and Mode Shapes 

The fundamental frequencies of the baseline NPSCuL FEM with conical isolators are 
shown in Table 5.6 along with the measured frequencies from the sine sweep test. The 
corresponding mode shapes are shown in Figure 5.4 and Figure 5.5. Similar to the base¬ 
line NPSCuL FEM, there is good correlation for the first two rocking modes, but the axial 
mode of the FEM is not as stiff as the measured frequency. However, this discrepancy was 
expected due to the poor correlation in the baseline NPSCuL FEM. Frequency correlation 
was not expected near and above 100 Hz, either, but it was not a concern because the higher 
modes do not typically affect the CLA results. 
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Table 5.6: Baseline NPSCuL with Conical Isolators — Fundamental Fre¬ 
quencies [1] 


Mode 

Frequencies (Hz) 

Percent 

Difference (%) 

FEM 

Test 

1 (Rocking, Y) 

20.3 

19.8 

2.42 

2 (Rocking, X) 

20.5 

20.4 

0.490 

3 (Axial, Z) 

42.7 

56.3 

24.2 

4 (Rocking/Bending, X) 

73.0 

77.8 

6.14 

5 (Rocking/Bending, Y) 

71.4 

81.5 

12.4 

6 (Shear/Rocking, Y) 

105 

144 

27.1 

7 (Shear/Rocking, X) 

103 

154 

33.1 
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(a) Baseline NPSCuL with Conical Isolators — 
Rocking Mode, X-Axis 


(b) Baseline NPSCuL with Conical Isolators 
Rocking Mode, Y-Axis 





(c) Baseline NPSCuL with Conical Isolators — Ax¬ 
ial "Pogo" Mode, Z-Axis 

Figure 5.4: Baseline NPSCuL with Conical Isolators — Fundamental Mode 
Shapes 
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(a) Baseline NPSCuL with Conical Isolators — 2nd (b) Baseline NPSCuL with Conical Isolators — 2nd 
Rocking Mode, X-Axis Rocking Mode, Y-Axis 




(c) Baseline NPSCuL with Conical Isolators — (d) Baseline NPSCuL with Conical Isolators — 
Shear Mode, X-Axis Shear Mode, Y-Axis 

Figure 5.5: Baseline NPSCuL with Conical Isolators — Secondary Mode 
Shapes 


5.3 Isogrid NPSCuL FEM with Conical Isolators 

The isogrid NPSCuL FEM, which consists of shell elements for the face sheets and beam 
elements for the stringers and fasteners [19], was modified to reflect the test conflguration 
with conical isolators; this model is shown in Figure 5.6. The baseplate, the flat ring, and 
the conical isolator portions of the FEM from the baseline NPSCuL FEM replaced the 
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unibase in the isogrid NPSCuL FEM. The fundamental frequencies of the isogrid NPSCuL 
FEM with conical isolators in each axis are shown in Table 5.7 along with the measured 
frequencies from the sine sweep test. 



Figure 5.6: Isogrid NPSCuL FEM with Conical Isolators 


Table 5.7: Isogrid NPSCuL with Conical Isolators — Fundamental Fre¬ 
quencies 


Mode 

Frequencies (Hz) 

Percent 

Difference (%) 

FEM 

Test 

1 (Rocking, Y) 

19.4 

19.5 

0.564 

2 (Rocking, X) 

19.6 

19.1 

2.56 

3 (Axial, Z) 

44.0 

58.4 

24.6 


The frequencies of the baseline and isogrid FEMs with conical isolators, shown in Ta¬ 
ble 5.8, confirm that the isogrid and baseline NPSCuL designs have the same fundamental 
modes in each axis (see Figure 4.6). The discrepancies for the higher modes are attributed 
to the stiffer walls and baseplate on the isogrid NPSCuL; however, the mode shapes are the 
same despite the differences in frequencies. 
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Table 5.8: Frequencies of Baseline and Isogrid NPSCuL FEMs with Conical 
Isolators 


Mode 

Frequencies (Hz) 

Percent 

Difference (%) 

Baseline 

FEM 

Isogrid 

FEM 

1 (Rocking, Y) 

20.3 

19.4 

4.43 

2 (Rocking, X) 

20.5 

19.6 

4.59 

3 (Axial, Z) 

42.7 

44.0 

3.0 

4 (Rocking/Bending, X) 

73.0 

83.8 

12.9 

5 (Rocking/Bending, Y) 

71.4 

86.5 

17.4 

6 (Shear/Rocking, Y) 

105 

133 

21.2 

7 (Shear/Rocking, X) 

103 

133 

22.9 


5.4 Baseline NPSCuL FEM with Cupmount Isolators 

The baseline NPSCuL FEM was also modified to reflect the test configuration with cup- 
mount isolators. Similar to the FEMs with the conical isolators, the baseline adapter ring 
was replaced with a flat ring and the isolators, which are modeled using thin-shell, con¬ 
centrated mass, beam, and spring elements (CBUSH elements in Nastran), respectively. 
The thinner sections of the ring that interface with the ABC plate are kept the same as the 
baseline adapter ring flange at 0.432 cm (0.17"), but the thicker sections of the ring that 
accommodate the isolators are 0.635 cm (0.25"); these thicknesses are modeled in the shell 
properties. The spring stiffness value in all degrees of freedom was 947 kN/m (5400 Ibf/in), 
as specified by the vendor’s isolator sizing analysis [37], which is performed by the vendor 
as a rough estimate of the isolator performance by modeling the isolated component as a 
single DOF at the CG of the component. The beam elements represent the base of each 
isolator, which does not contribute to the spring stiffness but provides a load path between 
the isolator and the adapter ring. The modulus of elasticity of these beam elements was ad¬ 
justed to 5.52 GPa (8.0x10^ psi) to represent the joint stiffness for correlation with the test 
frequencies. The non-structural mass of the isolators are accounted for in the concentrated 
mass elements located at the interface node between the beam and spring elements. The 
remaining NPSCuL structure and components were not modified. The FEM for the adapter 


105 




ring and isolators is shown in Figure 5.7. 



Figure 5.7: Flat Ring and Conical Isolator FEM 


5.4.1 Fundamental Modes — Frequencies 

The fundamental frequeneies of the baseline NPSCuL FEM with cupmount isolators are 
shown in Table 5.9 along with the measured frequeneies from the sine sweep test. Similar 
to the baseline NPSCuL FEMs with and without conical isolators, there is good correlation 
for the rocking modes, but the axial mode of the FEM is not as stiff as the measured 
frequency. The corresponding mode shapes are similar to those shown in Eigure 5.4 and 
Figure 5.5. 
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Table 5.9: Baseline NPSCuL with Cupmount Isolators — Fundamental 
Frequencies 


Mode 

Frequencies (Hz) 

Percent 

Difference (%) 

FEM 

Test 

1 (Rocking, Y) 

13.5 

12.4 

9.20 

2 (Rocking, X) 

13.7 

12.5 

8.93 

3 (Axial, Z) 

30.2 

39.5 

23.6 

4 (Rocking/Bending, X) 

53.9 

62.8 

14.2 

5 (Rocking/Bending, Y) 

54.3 

65.4 

17.0 

6 (Shear/Rocking, Y) 

80.3 

82.2 

2.28 

7 (Shear/Rocking, X) 

81.0 

85.6 

5.36 


5.5 Loads Analyses 

Loads analyses were only performed on the baseline NPSCuL with cupmount isolators 
because it is the most likely configuration to fly on a future ABC mission. In all axial load 
cases, the loads were applied to NPSCuL in tension (+Z direction). The structural integrity 
of the structure was also demonstrated during the qualification test. 

5.5.1 Fastener Analysis 

Similar to previous ABC missions, the static loads of 5 Gs lateral (X and Y) and 7 Gs axial 
(Z), as described in 3.1.1.4.2, were applied simultaneously to the FEM. The loads were 
extracted from the beam elements that were used to represent the fasteners. The maximum 
root-sum-square (RSS) shear and maximum axial loads at each joint type were used to cal¬ 
culate the margins of safety (MS); the NPSCuL fastener analysis tool was used to perform 
these calculations [19]. The margins of safety are shown in Table 5.10 and are positive at all 
joints, indicating that there will be no structural failures if the appropriate torque values are 
used. The torque values shown in Table 5.11, except for the ring/baseplate and ring/isolator 
joints, are calculated using the fastener analysis tool as well. The torque values for these 
two joints are instead based on experience with those bolt sizes and engineering judgment 
because the torques predicted by the analysis tool were higher than the typical range. The 
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values shown do not include running torque, but the final installation torques were used 
during the qualification test with no issues. 


Table 5.10: Margins of Safety, Fastener Analysis — Cupmount Isolators on 
NPSCuL 


Failure 

Joint 

Mode 

Ring/ 

ARC 

Baseplate/ 

Bracket/ 

PPOD/ 

SAD/ 

SAD 

Plate/ 

Ring/ 

Ring/ 


Plate 

Wall 

Wall 

Wall 

Wall 

SAD 

Frame 

Baseplate 

Isolator 

Yield 

0.10 

0.10 

0.11 

0.12 

0.11 

0.13 

0.75 

1.15 

Ultimate 

0.69 

0.72 

0.73 

0.73 

0.74 

0.76 

1.70 

2.29 

Shear 

3.2 

13 

4.4 

10 

18 

14.3 

7.2 

46.6 

Gapping 

1.20 

9.5 

4.4 

1.5 

21.5 

13.7 

1.9 

0.8 

Tear-Out 

1.9 

1.9 

0.52 

0.52 

0.53 

1.2 

0.1 

0.1 


Table 5.11: Torque Values For Baseline NPSCuL with Cupmount Isolators 


Joint 

Ring/ 

ABC 

Plate 

Baseplate/ 

Wall 

Bracket/ 

Wall 

PPOD/ 

Wall 

SAD/ 

Wall 

SAD 

Plate/ 

SAD 

Frame 

Ring/ 

Baseplate 

Ring/ 

Isolator 

Torque, 

in-lbs 

94 

101 

41 

39 

42 

12 

125 

50 


5.5.2 Stress Analysis 

Similar to the fastener loads analysis, the static loads requirement detailed in 3.1.1.4.2 was 
applied to the FEM. However, the adapter ring was modeled using solid brick elements for 
more accurate stress prediction. This model is shown in Figure 5.8. 
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Figure 5.8: Flat Ring FEM — Solid Brick Elements 


Factors of safety (FS) of 1.25 and 1.4 were applied to the resulting stresses; these are the 
most conservative factors amongst the common aerospace industry testing standards [21]. 
Yield and ultimate stress allowables of 393 kPa (57 ksi) and 469 kPa (68 ksi), respectively, 
were used to evaluate the margins [38]. The MS is calculated as follows: 


MS = 


^allowable 

(ypEM 


(5.1) 


where (Jaiiowabie is the allowable stress and Gpem is the maximum elemental Von Mises 
stress obtained from the FEM. As shown in Table 5.12, there are positive margins through¬ 
out the structure. The maximum stress is in the adapter ring at the ABC interface; this is 
shown in Figure 5.9. 

Table 5.12: Margins of Safety, Stress Analysis — Cupmount Isolators on 
NPSCuL 


Load Case 

Description 

Max Stress 

kPa [ksi] 

Yield 

MS 

Ultimate 

MS 

5 Gs Lateral (X,Y), 

7 Gs Axial (Z) 

182 [26.5] 

0.72 

0.83 
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1.83E+008 



1.67E+008 
1.52E+008 
1.37E+008 
1.22E+008 
1.07E+008 
9.14E+007 
7.62E+007 
6.10E+007 
4.58E+007 jjj 
3.06E+007 
1.54E+007 
2.11E+005 
Units = Pa(N/m''2) 


Figure 5.9: Von Mises Stress Distribution Due to ABC Static Loads 
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CHAPTER 6: 
Conclusion 


The combined use of FLVT and COTS isolators effectively reduces vibration levels for 
small satellites. This method, with the cupmount isolators on NPSCuL subjected to the 
ABC profile in particular, resulted in lowering the vibration Grms levels by 64 to 78% 
over the entire test frequency range. In the high-frequency range, the Grms responses were 
97% lower, and the magnitude of the P-POD acceleration response at the P-POD interface 
was 20 dB lower (2 orders of magnitude). Additionally, the cost of this modification at 
approximately $1,000 is about 2% of existing whole-spacecraft isolation systems. 

The incorporation of COTS isolators on to NPSCuL and using an FLVT set-up has also 
been demonstrated to be effective for a worst-case environment such as NPSCuL. It was 
shown that various COTS isolators can be used depending on the LV mounting configura¬ 
tion; the cupmount isolators are appropriate for any static loading condition, and the conical 
isolators are appropriate for mainly compressive static loading conditions. Table 6.1 sum¬ 
marizes how the design requirements have been met in modifying NPSCuL to incorporate 
cupmount isolators; this includes the completion of a successful qualification test. 

It was found that FLVT is more effective with vibration profiles that have high ASD values 
near the fundamental frequencies of the system, such as the ABC profile and other common 
vibration environments that are applicable to spacecraft that are similar in size to NPSCuL. 
It was also shown that COTS isolators are beneficial to small satellites, independent of a 
FLVT set-up. 
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Table 6.1: Launch Vehicle Requirements Summary 


Requirement 

Isolator/Re-Designed NPSCuL 

Feature 

AP static envelope must be within the 
auxiliary payload (AP) volume 

AP static envelope is the same in X and 

Y, 1.78 cm (0.7 in) reduction in Z 

AP mechanical interface must be 

compatible with the ABC plate 

Uses 8 of 24 available fasteners 

AP mass properties must be within the 
interface control document (ICD) range 

Increases NPSCuL mass by 1.13 kg (2.5 
lbs); no significant impact to AP mass 
properties 

AP must be able to withstand thermal 

environment 

Thermal environment is enveloped by 
operating range of elastomers 

AP must be able to withstand static loads 

environment 

Positive margins infastener and strength 
analysis 

AP must be able to withstand random 

vibration environment 

Baseline NPSCuL withcupmount 
isolators have been qualified to ABC 

levels 

AP fundamental frequency must be 
above 35 Hz or have no negative impacts 
to the LV or primary spacecraft 

Isolated NPSCuL fundamental 

frequenciesare below 35 Hz, but CLA 
results indicate no negative impacts to 

LV or primary spacecraft 


6.1 Summary of Contributions 

A novel, practical method of combining the use of COTS isolators and FLVT to reduce 
the vibration levels for small satellites has been developed. This method has not been 
previously used on small satellites, and it is compatible with the ESPA interface. The design 
can also be easily modified for different spacecraft interfaces and masses. The application 
of FLVT has also been extended to include the use of the semi-empirical method with 
isolated small spacecraft, which has now been validated. 

The reduction in vibration levels allows more sensitive and complex payloads to gain access 
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to space, which consequently improves the government’s infrastructure for launching small 
satellites because NPSCuL is the only government-owned launch vehicle interface that is 
also capable of launching multiple satellites from one auxiliary payload. Additionally, 
the reduction in cost by approximately 98% of existing whole-spacecraft isolation systems 
allows this method to be used on low-budget missions. 


6.2 Recommendations for Future Work 

Several topics may be investigated in the future to improve the isolator selection process 
and installation. An optimization program can be developed to determine the ideal isolator 
stiffness from a spacecraft’s mass, CG, and fundamental frequency requirements. Using 
cupmount isolators with non-threaded cores would improve the installation process by al¬ 
lowing a large torque value to be applied to the fastener without spinning the threaded 
core. Other types of secondary backout prevention would also be available, such as staking 
compound, instead of just a thread-lock adhesive. Capturing time histories during random 
vibration testing would also be beneficial for interpreting results. 

To further reduce the vibration test environment, the roll-off slope and value can be 
adjusted when flight data becomes available over the entire test frequency range. The 
data would allow for efforts to verify if a less conservative force limit would still provide 
adequate margin to actual flight levels. 
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APPENDIX A: 

Shock and Acoustic Testing No-Test Rationale 


Shock and acoustic environments were not driving requirements in re-designing NPSCuL, 
and no-test rationale was provided for the baseline NPSCuL configuration. As mentioned 
in [1], these environments and the analyses for the no-test rationales were performed to take 
the significant vibration reduction into account. It is noted that these will require formal 
approval from ULA in the event that the isolated NPSCuL is manifested for an auxiliary 
payload mission on the ABC. 

A.l Shock No-Test Rationale 

The ABC shock and vibration profiles are defined at the ABC-to-NPSCuL interface [2], 
the shock level is below the 50 in/sec velocity no-test line as specified in SMC-S-016 [12]. 
The shock profile is not intended to change for future missions; therefore, shock testing is 
still not required by ULA. 

A.2 Acoustic No-Test Rationale 

Similarly, acoustic analysis was performed assuming the ABC acoustic environment [2], 
shown in Figure A.l and defined at the ABC-to-NPSCuL interface, will not change for 
future missions. 
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Aft Bulkhead Carrier (ABC) Acoustic Requirement 



Frequency - Hz, 


I -' 3 OB Center 
Freqiienc^^ (Hz) 

Soiuid Pressure 
Level (dB) 

1.''3 OB C: enter 
Frequency (Hz) 

Sound Pressure 
Level (dB) 

32 
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40 
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50 
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63 
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117 5 
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127 
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117 
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129 

6300 

116 5 
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131 
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116 
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OiSPL = 

: 139.5 dB 



Figure A.l: Atlas V Maximum Predicted Acoustic Levels [2] 


The method described in a study performed by JPL [39] was used to determine that acous¬ 
tic testing was not required because the random vibration environment encompassed the 
acoustic levels; if the break-even surface area is higher than the component surface area, 
then the random vibration test will be more effective. Additionally, the study showed that 
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if the area/mass ratio of an object is less than 0.213 np-/kg (150 irp/lbm), the random 
vibration environment usually encompasses the acoustic profile. 

This analysis has been performed for all previous missions, and it was redone for 2 P-PODs 
or 1 6U dispenser mounted to the wall of NPSCuL using the isolated P-POD response level. 
The surface area calculated for both configurations was 0.152 (236 irp); this includes 

the conservative assumption that the entire NPSCuL wall is susceptible to acoustic loading. 
The total mass of the 6U configuration used in this analysis was 11.3 kg (25 Ibm), and the 
total mass of the two P-PODs was assumed to be 15.9 kg (35 lbs); both are conservatively 
low estimates. Although the sound pressure level (SPL) at 100 Hz, the estimated frequency 
of the panel mode, is 123.8 dB, the overall SPL (OASPL) of 139.5 dB was used for this 
analysis. The resulting break-even surface areas are 0.155 and 0.217 (240 and 337 irp) 

for the P-POD and 6U configurations, respectively. Although the break-even surface area 
is now lower for the isolated NPSCuL when compared to the baseline NPSCuL, it is still 
greater than the component surface area. Also, the break-even surface area/mass ratio of 
an NPSCuL wall with two P-PODs (using the lowest PSD value of 0.003 /Hz in the X- 
axis) is approximately 0.0137 nP/kg (9.6 irP/lb), which also indicates that acoustic testing 
is not required. The Matlab code used to perform this analysis is included in Appendix G. 
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APPENDIX B: 

P-POD Response Comparison Plots 




(a) X-Axis Envelope (b) Y-Axis Envelope 



(c) Z-Axis Envelope 

Figure B.l: Acceleration Response at P-POD Interface Due to ABC Profile, 
Configuration 1, with and without FLVT 
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(zh/,B) asv 



Config 2, FLVT — Config 2, Unnotched 

(a) X-Axis Envelope 



Config 2, FLVT—Config 2, Unnotched 

(b) Y-Axis Envelope 



(c) Z-Axis Envelope 

Figure B.2: Acceleration Response at P-POD Interface Due to ABC Profile, 
Configuration 2, with and without FLVT 
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ASD (g^/Hz) 



Config 3, FLVT—Config 3, Unnotched 



(a) X-Axis Envelope 


(b) Y-Axis Envelope 



(c) Z-Axis Envelope 

Figure B.3: Acceleration Response at P-POD Interface Due to ABC Profile, 
Configuration 3, with and without FLVT 
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APPENDIX C: 

Individual P-POD Responses, Grms 


This appendix contains all the Grms data for each P-POD by test axis for reference. Not 
all P-PODs were instrumented for every test configuration and input, but at a minimum, 
P-PODs 1 through 4 were instrumented to capture any trends due to differences in P-POD 
location. Only the data obtained for Configurations 1 through 3 due to the ABC input with 
and without FLVT are presented here. 


Table C.l: Configuration 1 P-POD Responses Due to ABC Input with 
FLVT 


P-POD 

Grms By Axis 

(Includes Cross-Axis Response) 

X 

Y 

Z 

1 

9.56 

6.51 

6.29 

2 

11.85 

8.12 

6.86 

3 

6.36 

10.79 

6.11 

4 

7.73 

12.58 

6.73 

5 

10.27 

6.14 

6.23 

6 

12.10 

7.36 

6.98 

7 

5.98 

7.75 

6.32 

8 

6.60 

10.53 

6.58 


123 





Table C.2: Configuration 1 P-POD Responses Due to ABC Input without 
FLVT 


P-POD 

Grms By Axis 

(Includes Cross-Axis Response) 

X 

Y 

Z 

1 

17.82 

9.94 

22.33 

2 

26.72 

11.87 

22.81 

3 

10.16 

17.50 

23.21 

4 

13.21 

28.20 

24.09 


Table C.3: Configuration 2 P-POD Responses Due to ABC Input with 
FLVT 


P-POD 

Grms By Axis 

(Includes Cross-Axis Response) 

X 

Y 

Z 

1 

3.11 

2.04 

3.72 

2 

4.46 

2.25 

3.78 

3 

2.24 

3.05 

3.79 

4 

2.34 

4.63 

3.91 

5 

3.38 

2.15 

3.90 

6 

5.05 

2.30 

3.83 
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Table C.4: Configuration 2 P-POD Responses Due to ABC Input without 
FLVT 


P-POD 

Grms By Axis 

(Includes Cross-Axis Response) 

X 

Y 

Z 

1 

5.15 

4.18 

5.71 

2 

7.66 

4.31 

5.69 

3 

3.64 

5.77 

6.68 

4 

3.66 

9.17 

6.93 

5 

5.32 

4.22 

5.95 

6 

8.06 

4.37 

5.84 


Table C.5: Configuration 3 P-POD Responses Due to ABC Input with 
FLVT 


P-POD 

Grms By Axis 

(Includes Cross-Axis Response) 

X 

Y 

Z 

1 

2.08 

1.18 

2.54 

2 

3.10 

1.29 

2.57 

3 

1.49 

1.86 

2.70 

4 

1.56 

2.72 

2.74 

5 

3.40 

1.48 

2.59 

6 

2.76 

1.60 

2.58 

7 

1.25 

2.20 

2.56 

8 

1.28 

3.11 

2.58 
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APPENDIX D: 

MATLAB Code — Force Limit Validation 


This appendix contains the MATLAB scripts for calculating the transmissibility of the base¬ 
line and isolated NPSCuL in a fixed-base and FLVT set-up; NPSCuL is modeled as a 
single degree-of-freedom in both cases. However, in the FLVT set-up, a second degree-of- 
freedom is introduced to explain the theory behind the notched input. 


%FLVT Illustration 

%Calculate accleration of mass and force in spring vs frequency 
%FLVT and standard random vibe set up 
%Using relative support motion equations 

o, 

o 

%W. Lan 

clc; clear all; close all; 

%SYSTEM PROPERTIES 
mass_sc = 200; %lbs 
mass_ABC = 20; %lbs 

k_LV = 2587221; 

% k_LV = 2737800; %axial K 

zheta = 0.02; 

%INPUT PARAMETERS 
omega = 0:2:3200; %Hz 

ABC_interface_MPE = [20 40 240 2000; 0.03 0.125 0.125 0.003]; 
ABC_MPE_full = 10.^interpl(logl0(ABC_interface_MPE(1^ :)),logl0( ... 
ABC_interface_MPE(2,:)),loglO(omega)); %ASD, g^2/Hz 

%CASE 1 - "EIXED BASE" CONEIGURATION 

%ASSUME SUPPORT MOTION USING ABSOLUTE GOORDINATES 
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%natural frequencies^ rad/s 
wn = sqrt(k_LV/mass_sc) ; 

disp( ’Fixed Base Model Damped Frequencies, Hz’) 
wdl = wn^sqrt(l-zheta^2) 

Y = 1; %sine sweep, g’s 
for i = 1:length(omega) 

%Displacement 

X(i) = Y^sqrt(1+(2^zheta^(omega(i)/wn))^2)/sqrt((1-(omega(i)/wn) . 

^2)^2+(2^zheta^omega(i)/wn)^2); 

%Acceleration 

% Xdoubledot (i) = X (i) Iomega (i) ^^2/386.4; 

end 

%Acceleration of mass due to 7.6 Grms input 

accel = X.^ABC_MPE_full; 

%Force in Spring 
k_LV_force = accel^mass_sc; 

%CASE 2 - FLVT CONFIGURATION 

%ASSUME VIBRATION ABSORBER MODEL%ABC PLATE IS THE ABSORBER (CONTRARY 
TO BC 
%SETUP) 

%"2" refers to NPSCuL 
%"1" refers to ABC plate 
ml = mass_ABC; 
m2 = mass_sc; 

kl = k_LV; %assume fixed base stiffness 
k2 = 2737800; %axial K 

disp(’FLVT Model Damped frequencies, Hz’) 
wl = sqrt(kl/ml) 
w2 = sqrt(k2/m2) 

mu = m2/ml; 
q = w2/wl; 

toggle =2; %1 for undamped, 2 for damped 
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if toggle 


1 


%FROM KELLY, EQN 8.42 AND 8.43 - UNDAMPED VIBRATION ABSORBER 
for i = 1:length(omega) 
rl = omega(i)/wl; 
r2 = omega(i)/w2; 

Ami(i) = abs((l-r2^2)/(rl^2^r2^2-r2^2-(1+mu)^rl^2+l)); 
Am2(i) = abs (1/(rl^2^r2^2-r2^2-(1+mu)^rl^2 + l)); 
if Ami(i) == 0 

Aml(i) = lE-4; 

end 

if Am2(i) == 0 

Am2(i) = lE-4; 

end 

end 

else 


%EROM KELLY, EQN 8.54 AND 8.55 - DAMPED VIBRATION ABSORBER 
for i = 1:length(omega) 
rl = omega(i)/wl; 
r2 = omega(i)/w2; 

Aml(i) = sqrt ( ( (2^zheta^rl^q) ^2+(rl^2-q^2) "^2) / ( (rl^4-(1+(1+mu 
) ^q^2) ^rl^2+q'^2) ^2 + (2^zheta^rl^q) ^2^ (l-rl^2^ (1+mu) ) ^2) ) ; 
Am2(i) = sqrt((q^4+(2^zheta^q)^2)/((rl^4-(1+(1+mu)^q^2)^rl^2+ 
q^2) ^2+(2^zheta^rl^q)^2^(l-rl^2^(1+mu) )^2) ) ; 
if Ami(i) == 0 

Aml(i) = lE-4; 

end 

if Am2(i) == 0 

Am2(i) = lE-4; 

end 

end 


end 


%Acceleration of S/C 
accel2 = Am2.^ABC_MPE_full; 
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%Force in Spring^ S/C to LV 
kforceS = accel2^mass_sc; 

%Acceleration of ABC plate 
accell = Ami.^ABC_MPE_full; 

%Force in Spring, S/C to LV 
kforcel = accell^mass_ABC; 

%TRANSMISSIBILITY PLOT 
figure (’ Color' , [1 1 1]) 
loglog(omega,X) 
hold on 

xlim([20 2000]) 
grid on 

ylabel( ’Transmissibility' , 'FontSize' ,24) 

loglog(omega,Ami) 

plot ([117,117], [lE-4,1E3] , 'k—' ) 

xlim([20 2000]) 

grid on 

xlabel (’ Frequency, Hz', 'FontSize' ,24) 

legendl=legend( 'SDOF,NPSCuL' , '2DOF,ABC Plate (m_l)' , 'Location' , ' ... 

SouthOutside' , 'Orientation' , 'Horizontal' ); 
set(legendl, 'FontSize' ,24) 

% title('Vibration Absorber Effect in a 2DOF Model' , 'FontSize' , 24) 

%2DOF (FLVT) PLOT - % Acceleration of ABC plate and test control 
load testdata.mat 
figure (' Color' , [1 1 1]) 

loglog(omega,accell, '-b' , 'LineWidth' ,1.5) 
hold on 

loglog(freq,configl.ABC.control.Z, '-g', 'LineWidth ', 1.5) 
xlim([20 2000]) 
grid on 

xlabel (' Frequency, Hz', 'FontSize' ,24) 

% title('2DOF Model Comparison with Force Limited Test Control','... 
FontSize',24) 

ylabel (' Acceleration Spectral Density (g^2/Hz)' , 'FontSize' ,24) 

% legend('Fixed Base (NPSCuL)','FLVT,ml(ABC Plate)','FLVT,m2(NPSCuL)... 
','Location','SouthOutside','Orientation','Horizontal') 


130 


Notched ... 


Iegend2=legend (' Acceleration Response of ABC Plate' 

Acceleration Input (Test)'Location','SouthOutside' , 'Orientation. 
' , 'Horizontal' ); 
set(legendC, 'FontSize' ,24) 

%SiMPLE TDFS METHOD 
Q=l/(2^zheta) ; 
c=2^zheta; 

betal = 1+1/2^mu+sqrt(mu+1/4^mu^2); 
beta2 = l+l/2^mu-sqrt(mu+l/4^mu^2); 
for i = 1:251 

% H2(i)=-(k2 + c^omega (1) ^ j )/(kl ^m2 Iomega (i)'^2 - kl^k2 + k2^ml^.. 

omega (1)^2 + k2^m2Iomega (1)'^2 - ml^m2Iomega (1) 4 - c^klIomega (1) ^ . 

j + c^ml^omega (1) "^3^ j + c^m2^omega (1) ^^3^ j ) ; 

% Sff(i) = (l + (omega (1) /w2/Q) ^^ 2 ) / ( (1- (omega (1) /w2) ^^2) ^2 + (omega (1) / . 

w2/Q)^2)^ABC_MPE_full(1)^abs(H2(1))^2; 

Sff(1) = ABC_MPE_full(1)^m2^2^(1+(omega(1)/w2)^2/Q^2)/((1-(omega( 
1)/w2)^2)^2+(omega(1)/w2)^2/Q^2)^4.448221628254617; 

"o 

end 

Sffl = (l+betal/Q^2)/((1-betal)^2+betal/Q^2)^ABC_MPE_full^m2^2; 

Sff2 = (l+beta2/Q^2)/((l-beta2)^2+beta2/Q^2)^ABC_MPE_full(1:251)^m2 .. 
-2; 

Sff3 = [1.52E+03,2.60E+03,6.31E+03,5.57E+03,2.95E+03,2.95E+03,2.59E .. 
+03,1.15E+03,7.06E+01]^4.448221628254617; 

%COMPLEX TDES METHOD 
%from Excel spreadsheet 

f_center = [20,25,40,64,102,162,258,410,500]; 

Sff_complex = [1.07E+03,1.84E+03,4.46E+03,1.58E+04,1.48E+04,4.46E ... 
+03,3.93E+03,1.74E+03,1.07E+02]^4.448221628254617; 

%MAX OF BOTH 

Sff_max = max([Sff_complex;Sff3]); 

%SEM1-EMP1R1CAL METHOD 

[E,T,dataout] = forcelimitcalc(ABC_interface_MPE',w2,500,0,m2); 
%C0MPAR1S0N PLOT, FORCE LIMIT METHODS 
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figure (’ Color ',[ 1 1 1]) 

% loglog(omega(1:251)^abs(Sff2)^ omega(1:251)^abs(Sff)^ dataout ( :^1)^ . 

dataout(:,2),f_center,Sff_complex^ f_center,Sff3) 

% loglog(dataout {:,!), dataout(:,2),f_center^ Sff_complex,f_center,Sff3 

) 

% loglog(dataout(:^1)^ dataout(:^ 2)^ f_center^ Sff_max) 

loglog(dataout(:,1)^ dataout(:^2)^4.448221628254617^ 'k-', 'LineWidth' . . 
.2) 

hold on 

scatter(f_center,Sff_complex, 40, ’r’, 'filled’) 
scatter(f_center,Sff3,40, 'b', 'filled') 
loglog(f_center,Sff_max, ’g-' , 'LineWidth ', 2) 
xlim([20 2000]) 
grid on 

xlabel (’ Frequency, Hz’, ’FontSize’ ,24) 

ylabel (’ Force Spectral Density (N^2/Hz) ’ , ’FontSize’ ,24) 

% ylabel ('Force Spectral Density (lb^2/Hz) ', 'FontSize',24) 
set(gca, ’FontSize' ,24) 

% legend2=legend('Semi-Empirical', 'TDFS', 'Location', 'SouthOutside' , ' . 
Orientation','Horizontal'); 

legend2=legend (' Semi-EmpiricalComplex TDES ',' Simple TDES',’TDES .. 

Max' , 'Location' , 'SouthOutside’ , ’Orientation’ , 'Horizontal' ); 
set(legends, 'EontSize' ,24) 


%ELVT Illustration, Isolated NPSCuL 

%Calculate accleration of mass and force in spring vs frequency 
%ELVT and standard random vibe set up 
%Using relative support motion equations 

'o 

%W. Lan 

clc; clear all; close all; 

%SYSTEM PROPERTIES 
mass_sc = 189; %lbs 
mass_ABC = 20; %lbs 
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% k_LV = 605000; %lb/in, lateral 
k_LV = 2737800; %axial K 

zheta = 0.02; 

%1NPUT PARAMETERS 
omega = 0:1:2000; %Hz 

ABC_interface_MPE = [20 40 240 2000; 0.03 0.125 0.125 0.003]; 
%"ABC_MPE_full" is really just 0.5 G sine sweep to compare with ... 

apparent 
%mass #'s 

% ABC_MPE_full = 0.5; 

ABC_MPE_full = 10.^interpl(logl0(ABC_interface_MPE(1,:)),logl0( ... 

ABC_interface_MPE(2,:)),loglO(omega)); %ASD, g^2/Hz 
%plot ABC_MPE_full to check 
% figure 

% loglog(omega^ (ABC_MPE_full)) 

% xlim([20 2000]) 

% grid on 

%CASE 1 - EIXED BASE CONEIGURATION - NO ISOLATORS 
%ASSUME SUPPORT MOTION USING ABSOLUTE GOORDINATES 

%natural frequencies^ rad/s 
wn = sqrt(k_LV/mass_sc); 

disp('Eixed Base Model Damped frequencies^ Hz’) 
wdl = wn^sqrt(l-zheta^2) 

Y = 1; %sine sweep^ g's 
for 1 = 1:length(omega) 

%Displacement 

X(1) = Y^sqrt(1+(2^zheta^(omega(1)/wn))^2)/sqrt((1-(omega(1)/wn) ... 

^2)^2+(2^zheta^omega(1)/wn)^2); 

%Acceleration 

% Xdoubledot (1) = X (1) Iomega (1)'^2/386.4; 

end 

%Acceleration of mass due to 7.6 Grms input 
accel = X.^ABC_MPE_full; 
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%Force in Spring 
% k_LV_force = accel^mass_sc; 

%CASE 2 - FLVT CONFIGURATION - NO ISOLATORS 

%ASSUME VIBRATION ABSORBER MODEL 

%USING EQN 3.15 FROM HARTOG, P. 89 

%"2" refers to NPSCuL 

% II 111 2 -efers to ABC plate 

ml = mass_ABC; 

m2 = mass_sc; 

kl = k_LV; %assume fixed base stiffness 
% k2 = 460800; %lateral K 
k2 = 2737800; %axial K 

disp('FLVT Model Damped Frequencies, Hz’) 
wl = sqrt(kl/ml) 
w2 = sqrt(k2/m2) 


mu = m2/ml; 
q = w2/wl; 

toggle =2; %1 for undamped, 2 for damped 
if toggle == 1 

%FROM KELLY, EQN 8.42 AND 8.43 - UNDAMPED VIBRATION ABSORBER 
for i = 1:length(omega) 
rl = omega(i)/wl; 
r2 = omega(i)/w2; 

Aml(i) = abs((l-r2^2)/(rl^2^r2^2-r2^2-(1+mu)^rl^2+l)); 
Am2(i) = abs(1/(rl^2^r2^2-r2^2-(1+mu)^rl^2+l)); 
if Ami(i) == 0 

Aml(i) = lE-4; 

end 

if Am2(i) == 0 

Am2(i) = lE-4; 

end 

end 
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else 


%FROM KELLY, EQN 8.54 AND 8.55 - DAMPED VIBRATION ABSORBER 
for I = 1:length(omega) 
rl = omega(I)/wl; 
r2 = omega(I)/w2; 

Ami(I) = sqrt(((2^zheta^rl^q)^2+(rl^2-q^2)^2)/((rl^4-(1+(1+mu 
)^q^2)^rl^2+q^2)^2+(2^zheta^rl^q)^2^(l-rl^2^(1+mu))^2)); 
Am2 (1) = sqrt ( (q^4+(2^zheta^q) ^2 ) / ( (rl^4-(1+(1+mu)’^q^2 ) ^rl^2 + 

q^2) ^2+(2^zheta^rl^q)^2^(l-rl^2^(1+mu) )^2) ) ; 
if Ami(1) == 0 

Aml(i) = lE-4; 

end 

if Am2(1) == 0 

Am2(i) = lE-4; 

end 

end 

end 

%CASE 3 - ELVT CONEIGURATION - WITH ISOLATORS 
% k3 = 302400; %lateral K 

k3 = 2.9489e+05; %axial K to match sine sweep data 
k3 = 472500; %axial K to match random test data 
zheta = 0.22; 

%k3 = 99981; %isolator stiffness 

% wl = sqrt(kl/ml) 
w3 = sqrt(k3/m2) 

disp('FLVT Isolated Model Damped Frequencies, Hz') 
wl 

wd3 = w3^sqrt(l-zheta^2) 
q3 = w3/wl; 
if toggle == 1 

%FROM KELLY, EQN 8.42 AND 8.43 - UNDAMPED VIBRATION ABSORBER 
for 1 = 1:length(omega) 
rl = omega(i)/wl; 


135 


r3 = omega(i)/w3; 

Aml3(i) = abs ( (l-r3^2 ) / (rl^2^r3'^2-r3^2-(1+mu) ^rl^2 + 1) ) ; 
Am23(i) = abs(1/(rl^2^r3^2-r3^2-(1+mu)^rl^2+l)); 
if Ami3(i) == 0 

Aml3(i) = lE-4; 

end 

if Am23(i) == 0 

Am23(i) = lE-4; 

end 

end 


else 

%EROM KELLY, EQN 8.54 AND 8.55 - DAMPED VIBRATION ABSORBER 
for i = 1:length(omega) 
rl = omega(i)/wl; 

Aml3(i) = sqrt( ( (2^zheta^r1^q3)^2+(rl^2-q3^2) ^2) /( (rl . . . 

^4-(1+(1+mu)^q3^2)^rl^2+q3^2)^2+(2^zheta^r1^q3)^2^(1-rl ... 
"^2^ (1+mu) ) ^ 2 ) ) ; 

Am23(i) = sqrt((q3^4+(2^zheta^q3)^2)/((rl^4-(1+(1+mu)^q3^2)^ . 

rl^2+q3^2)^2+(2^zheta^rl^q3)^2^(l-rl^2^(1+mu))^2)); 
if Ami3(i) == 0 

Aml3(i) = lE-4; 

end 

if Am23(i) == 0 

Am23(i) = lE-4; 

end 

end 

end 

%EIXED BASE MODEL - WITH ISOLATORS 
Y = 1; %sine sweep, g's 
for i = 1:length(omega) 

%Displacement -if Y = 1, then same as amplitude ratio, X/Y 
X3(i) = Y^sqrt(1+(2^zheta^(omega(i)/w3))^2)/sqrt((1-(omega(i)/w3) 
^2)^2+(2^zheta^omega(i)/w3)^2); 

%Acceleration 

% Xdoubledot (i) = X (i) Iomega (i) "^2/386.4; 

end 
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%2D0F (FLVT) PLOT - % Acceleration of ABC plate and test control 
load testdata.mat 
figure (’ Color’ , [1 1 1]) 

loglog(omega,Aml3.^ABC_MPE_full, ’-b’ , ’LineWidth’ ,1.5) 
hold on 

loglog(freq,config2d.ABC.control.Z, '-g', 'LineWidth' ,1.5) 
xlim([20 2000]) 
grid on 

xlabel (’ Frequency, Hz’, 'FontSize’ ,24) 

% title(’2DOF Model Comparison with Force Limited Test Control, ... 
Isolated NPSCuL','FontSize',24) 

ylabel (’ Acceleration Spectral Density (g^2/Hz)' , 'FontSize' ,24) 

% legend(’Fixed Base (NPSCuL)’,’FLVT,ml(ABC Plate)’,’FLVT,m2(NPSCuL). 

’,’Location’,’SouthCutside’,’Crientation’,’Horizontal’) 
legend2=legend (’ Acceleration Response of ABC Plate ’,’ Notched ... 

Acceleration input (Test) ’, ’Location’,’SouthCutside’,’Crientation 
’ , ’Horizontal’ ); 
set(legends, ’FontSize ’, 24) 

%TRANSMiSSiBiLiTY PLCT 

%M1, FLVT PLATE 

figure (’ Color ’,[ 1 1 1]) 

loglog(omega,X3, ’-k’ , ’LineWidth’ ,1.5) 

hold on 

loglog(omega,Aml3, ’-b’ , ’LineWidth ’, 1.5) 
xlim([5 2000]) 
grid on 

xlabel (’ Frequency, Hz’, ’FontSize’ ,24) 
ylabel( ’Transmissibility’ , ’FontSize’ ,24) 

% title(’Vibration Absorber Effect in a 2DCF Model - isolated NPSCuL. 

’) 

legendl=legend( ’Fixed Base (SDCF), isolated NPSCuL ’,’ 2DCF ABC Plate , 
ml (isolated NPSCuL)’ , ’Location’ , ’SouthCutside’ , ’Crientation’ , ’ .. 
Horizontal’ ) 

set(legendl, ’FontSize ’, 24) 

% bigtitle = sprintf(’Effect of FLVT (2DCF Damped TMD Model) on ... 
Baseline and isolated NPSCuL Fundamental Frequencies wrt ABC ... 
Plate’); 
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suptitle(bigtitle) 
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APPENDIX E: 

MATLAB Code — Force Limit Calculation 


This appendix contains the MATLAB script and associated functions for calculating the 
force limit using the semi-empirical method. The inputs are the unnotched ASD values, the 
roll-off frequency, and the mass of the test article. The force limit is plotted and exported 
as a text file. 


clc; clear all; close all; 

%calculate force limit using the semi-empirical method 
load randomvibeinputs 
%input at MPE 

MPE = input (’Enter desired input name (see available inputs in ... 
workspace): ’ ); 

%scaled MPE 

scale = input (’Enter desired scale of MPE (dB): ’); 

MPE(:,2) = 10.^ (scale./lO) .^MPE(:,2) ; 

%roll-off frequency and ASD value 

freq_rolloff = input (’Enter roll-off frequency (Hz): ’); 

ASD_rolloff = logslope(MPE,freq_rolloff); 

%cut-off frequency and ASD value 
freq_cutoff = 500; 

ASD_cutoff = logslope(MPE,freq_cutoff); 

%new sorted MPE data 

ASD_forcelimit_all = sortrows([MPE;[freq_rolloff,ASD_rolloff];[ ... 
freq_cutoff,ASD_cutoff]]); 

ASD_forcelimit = ASD_forcelimit_all((find(ASD_forcelimit_all(:,!)<. .. 
freq_cutoff)),:); 
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%C value 

C = sqrt(2); 

%Mass of test article 

MO = input (’Enter mass of test article(lbs): '); 

[F,T,dataout] = forcelimitcalc(MPE,freq_rolloff,freq_cutoff^scale,MO) . 

} 

format longg 

csvwrite( ’Force Limit.dat’ ,dataout) 

for 1 = 1:length(ASD_forcelimit)-1 
y2 = dataout(1 + 1,2) ; 
yl = dataout(1, 2) ; 
x2 = dataout(1 + 1, 1) ; 
xl = dataout(1, 1) ; 
dB = lO^loglO(y2/yl); 

Oct = log(x2/xl)/log(2) ; 
slope (1) = dB/Oct; 

end 


function y = logslope(MPE,freq_rolloff); 

%interpolate y value of log scale 

%find slope at roll-off frequency 
%find xl,yl 

temp = MPE((find(MPE(:,1)<freq_rolloff)),:); 
[dump,idxl] = min(abs(freq_rolloff-temp(:,1))); 
xl = temp(idxl,1); 
yl = temp(idxl,2); 

%find x2,y2 

temp = MPE((find(MPE(:,1)>freq_rolloff)),:); 
[dump,idx2] = min(abs(freq_rolloff-temp(:,1))); 
x2 = temp(idx2,1); 
y2 = temp(idx2,2); 

%calculate slope 
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dB = lO^loglO(y2/yl); 

Oct = log(x2/xl)/log(2) ; 
slope = dB/Oct; 

%ASD value at roll-off frequency 

y = yl^(freq_rolloff/xl)^(slope/(lO^loglO(2))); 


function [F,T,dataout] = forcelimitcalc(MPE,freq_rolloff,freq_cutoff, 
scale,MO); 

%function to calculate force limit using the semi-empirical method 

%scaled MPE 

% MPE(:,2) = 10.^(scale./lO).^MPE(:,2); 

%roll-off frequency ASD value 
ASD_rolloff = logslope(MPE,freq_rolloff); 

%cut-off frequency ASD value 
ASD_cutoff = logslope(MPE,freq_cutoff); 

%new sorted MPE data 

ASD_forcelimit_all = sortrows([MPE;[freq_rolloff,ASD_rolloff];[ ... 
freq_cutoff,ASD_cutoff]]); 

ASD_forcelimit = ASD_forcelimit_all((find(ASD_forcelimit_all(:,!)<. .. 
freq_cutoff)) , :) ; 

%C value 
C = sqrt(2); 

%calculate force limit 

for 1 = 1:length(ASD_forcelimit); 

if ASD_forcelimit(1,1) < freq_rolloff 
force_limit(1,1) = C^2^M0^2^ASD_forcelimit(1,2); 
else % ASD_forcelimit(1,1) >freq_rolloff 

force_limit(1,1) = 0^2^M0^2^ASD_forcelimit(1,2)^(freq_rolloff 
/ASD_forcelimit(1,1))^2; 

end 

end 
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freq = ASD_forcelimit(:^1); 

F = figure; 

% plotyy(ASD_forcelimit(:^1),force_limit , ASD_forcelimit 
ASD_forcelimit(:,2)^ 'loglog','plot'); 
loglog(freq,force_limit, 'ks-') 
grid on 
hold on 

loglog(ASD_forcelimit_all(:,1),ASD_forcelimit_all(:,2), 'bs-') 
xlabel (' Frequency^ Hz') 

axis([20 2000 0.l^min(ASD_forcelimit(:^2)) 1.1^max(force_limit)]) 


format short 

T = table(freq,force_limit) 

dataout = [freq,force_limit]; 

% dataout = [freq;force_limit]; 
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APPENDIX F: 

MATLAB Code — Data Processing 


This appendix contains the MATLAB scripts and associated functions for reading in data 
from an Excel spreadsheet, finding the max ASD values for a specified P-POD position or 
test axis, and calculating the Grms over a specified frequency range. This appendix also 
contains the MATLAB scripts for generating figures from the measured data. 


clc; clear all; close all; 
global freq a! 

%Run vibedataread and plot data 

%cutoff frequency for low/high freq content^ Hz 
cutoff = 500; 

ABCfilename = 'NPSCuL_8P2M2s_ABCMPE+0dB' ; 

ABCUfilename = '131107_NPSCuL_8P2M2s_ABCMPE+0dB_Unnotched' ; 
ORSSfilename = 'NPSCuL_8P2M2s_ORS3MPE-6dB' ; 

ORS3Ufilename = 'NPSCuL_8P2M2s_ORS3MPE-6dB_Unnotched' ; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ABC MPE ... 

ooooooooooooooooooooooooooooooooooooooooooooooooooooo 

[freq,Af,config.ABC.ppod.X,config.ABC.force.X,config.ABC.control.X] = ... 
vibedataread(ABCfilename, 'X+OdB’ ); 

[- 1 ,-I, config. ABC . ppod. Y, config . ABC . force . Y, config . ABC . control. Y] = ... 
vibedataread(ABCfilename, ’Y+OdB’ ); 

[- 1 ,-I, config. ABC . ppod. Z, config . ABC . force . Z, config . ABC . control. Z ] = ... 
vibedataread(ABCfilename, 'Z+OdB' ); 


%calculate Grms for each P-POD in every axis for every test axis 
config.ABC.grms.X.all = grmscalc(20,2000, 'p’ ,config.ABC.ppod.X,8); 
config.ABC.grms.X.low = grmscalc(20,cutoff, 'p' ,config.ABC.ppod.X,8); 
config.ABC.grms.X.high = config.ABC.grms.X.all-config.ABC.grms.X.low; 
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config.ABC.grms.Y.all = grmscalc(20,2000, 'p' ,config.ABC.ppod.Y,8); 
config.ABC.grms.Y.low = grmscalc(20,cutoff, 'p' ,config.ABC.ppod.Y,8); 
config.ABC.grms.Y.high = config.ABC.grms.Y.all-config.ABC.grms.Y.low; 

config.ABC.grms.Z.all = grmscalc(20,2000, 'p’ ,config.ABC.ppod.Z,8); 
config.ABC.grms.Z.low = grmscalc(20,cutoff, 'p’ ,config.ABC.ppod.Z,8); 
config.ABC.grms.Z.high = config.ABC.grms.Z.all-config.ABC.grms.Z.low; 

%calculate Grms for envelope of each P-POD 

ppodnum = size(config.ABC.ppod.X.x,2); % # of P-PODs in data set - .. 

should be 8 unless it is unnotched ABC of baseline (old data) 
config.ABC.max.PPOD = maxASD(config.ABC.ppod,ppodnum); % max ASD ... 
values of each P-POD in each axis 

config.ABC.grmsall.all = grmscalc(20,2000, 'p', config.ABC.max.PPOD, 8); 
config.ABC.grmsall.low = grmscalc(20,cutoff, 'p' ,config.ABC.max.PPOD .. 
.8) ; 

config.ABC.grmsall.high = config.ABC.grmsall.all- ... 
config.ABC.grmsall.low; 

%calculate Grms for envelope of each axis 

config.ABC.max.Axis = maxAxis(config.ABC.max.PPOD); %max ASD values . 
of all P-PODs in each axis 

config.ABC.grmsAxis.all = grmscalc(20,2000, 'd' ,config.ABC.max.Axis,8) 

r 

config.ABC.grmsAxis.low = grmscalc(20,cutoff, ’d ', config.ABC.max.Axis . 

,8) ; 

config.ABC.grmsAxis.high = config.ABC.grmsAxis.all- ... 
config.ABC.grmsAxis.low; 


0000000000000000000000000000000000000000000000000000 wr\0 O JUj O Ullj • • • 

UNNOTCHED %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


[- 1 ,- 1 , conf ig. ORSSUnnotched. ppod. X, conf ig . ORSSUnnotched. force . X, . . . 

config.0RS3Unnotched.control.X] = vibedataread(ORSSUfilename, ’X-6 
dB'); 

[- 1 ,-I, conf ig. ORSSUnnotched. ppod. Y, conf ig . ORSSUnnotched. force . Y, . . . 

config.0RS3Unnotched.control.Y] = vibedataread(0RS3Ufilename, ’Y-6 
dB'); 

[- 1 ,- 1 , conf ig. 0RS3Unnotched. ppod. Z, conf ig . 0RS3Unnotched. force . Z, . . . 

config.0RS3Unnotched.control.Z] = vibedataread(0RS3Ufilename, 'Z-6 
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dB ' ) ; 


%calculate Grms for each P-POD in every axis for every test axis 
config.ORSSUnnotched.grms.X.all = grmscalc(20,2000, ’p’ , ... 

config.ORS3Unnotched.ppod.X,8); 
config.ORSSUnnotched.grms.X.low = grmscalc(20,cutoff, 'p . 
config.ORS3Unnotched.ppod.X,8); 

config.ORS3Unnotched.grms.X.high = config.ORS3Unnotched.grms.X.all- . 
config.ORS3Unnotched.grms.X.low; 

config.ORS3Unnotched.grms.Y.all = grmscalc(20,2000, ’p . 

config.ORS3Unnotched.ppod.Y,8); 
config.ORS3Unnotched.grms.Y.low = grmscalc(20,cutoff, ’p . 
config.ORS3Unnotched.ppod.Y,8); 

config.ORS3Unnotched.grms.Y.high = config.ORS3Unnotched.grms.Y.all- . 
config.ORS3Unnotched.grms.Y.low; 

config.ORS3Unnotched.grms.Z.all = grmscalc(20,2000, ’p . 

config.ORS3Unnotched.ppod.Z,8); 
config.ORS3Unnotched.grms.Z.low = grmscalc(20,cutoff, ’p . 
config.ORS3Unnotched.ppod.Z,8); 

config.ORS3Unnotched.grms.Z.high = config.ORS3Unnotched.grms.Z.all- . 
config.ORS3Unnotched.grms.Z.low; 

%calculate Grms for envelope of each P-POD 

ppodnum = size(config.ORS3Unnotched.ppod.X.x,2); % # of P-PODs in .. 
data set - should be 8 unless it is unnotched ABC of baseline ( . . 
old data) 

config.ORS3Unnotched.max.PPOD = maxASD(config.ORS3Unnotched.ppod, ... 

ppodnum); % max ASD values of each P-POD in each axis 

config.ORS3Unnotched.grmsall.all = grmscalc(20,2000, ’p ’,.. . 

config.ORS3Unnotched.max.PPOD,8); 
config.ORS3Unnotched.grmsall.low = grmscalc(20,cutoff ,' p ',.. . 
config.ORS3Unnotched.max.PPOD,8); 

config.ORS3Unnotched.grmsall.high = config.ORS3Unnotched.grmsall.all 
config.ORS3Unnotched.grmsall.low; 

%calculate Grms for envelope of each axis 
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config.ORSSUnnotched.max.Axis = maxAxis(config.ORSSUnnotched.max.PPOD 
); %max ASD values of all P-PODs in each axis 
config.ORSSUnnotched.grmsAxis.all = grmscalc(20,2000, ’d' , ... 

config.ORS3Unnotched.max.Axis,8); 
config.ORS3Unnotched.grmsAxis.low = grmscalc(20,cutoff ,' d . 

config.ORS3Unnotched.max.Axis,8); 
config.ORS3Unnotched.grmsAxis.high = ... 
config.ORS3Unnotched.grmsAxis.all- ... 
config.ORS3Unnotched.grmsAxis.low; 




ORS3 MPE-6dB ... 


ooooooooooooooooooooooooooooooooooooooooooooooooooooo 

[- 1 ,- 1 , config. ORS3 . ppod .X, config . ORS3 . force.X, config. ORS3 . control .X] = 
vibedataread(ORS3filename, ’X-6dB’ ); 

[- 1 ,- 1 , config. ORS3 . ppod. Y, config . ORS3 .force.Y,config. ORS3 . control. Y] = 
vibedataread(ORS3filename, ’Y-6dB’ ); 

[- 1 ,- 1 , config. ORS3 . ppod. Z, conf ig . ORS3 .force.Z,con fig. ORS3 . control. Z ] = 
vibedataread(ORS3filename, 'Z-6dB’ ); 


%calculate Grms for each P-POD in every axis for every test axis 
config.ORS3.grms.X.all = grmscalc(20,2000, ’p' ,config.ORS3.ppod.X, 8) ; 
config.ORS3.grms.X.low = grmscalc(20,cutoff, 'p' ,config.ORS3.ppod.X,8) 

r 

config.ORS3.grms.X.high = config.ORS3.grms.X.all- ... 
config.ORS3.grms.X.low; 

config.ORS3.grms.Y.all = grmscalc(20,2000, ’p’ ,config.ORS3.ppod.Y,8); 
config.ORS3.grms.Y.low = grmscalc(20,cutoff, 'p’ ,config.ORS3.ppod.Y,8) 

r 

config.ORS3.grms.Y.high = config.ORS3.grms.Y.all- ... 
config.ORS3.grms.Y.low; 

config.ORS3.grms.Z.all = grmscalc(20,2000, ’p' ,config.ORS3.ppod.Z,8); 
config.ORS3.grms.Z.low = grmscalc(20,cutoff, ’p’ ,config.ORS3.ppod.Z,8) 

r 

config.ORS3.grms.Z.high = config.ORS3.grms.Z.all- ... 
config.ORS3.grms.Z.low; 

%calculate Grms for envelope of each P-POD 
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ppodnum = size(config.ORSS.ppod.X.x,2); % # of P-PODs in data set - . 

should be 8 unless it is unnotched ABC of baseline (old data) 
config.0RS3.max.PPOD = maxASD(config.0RS3.ppod,ppodnum); % max ASD .. 

values of each P-POD in each axis 

config.0RS3.grmsall.all = grmscalc(20,2000, ’p’, config.ORS3.max.PPOD .. 

.8) ; 

config.ORS3.grmsall.low = grmscalc(20,cutoff, 'p ’, config.ORS3.max.PPOD 

.8) ; 

config.ORS3.grmsall.high = config.ORS3.grmsall.all- ... 
config.ORS3.grmsall.low; 


%calculate Grms for envelope of each axis 

config.ORS3.max.Axis = maxAxis(config.ORS3.max.PPOD); %max ASD values 
of all P-PODs in each axis 

config.ORS3.grmsAxis.all = grmscalc(20,2000, 'd ', config.ORS3.max.Axis . 

.8) ; 

config.ORS3.grmsAxis.low = grmscalc(20, cutoff ,' d ',.. . 
config.ORS3.max.Axis,8); 

config.ORS3.grmsAxis.high = config.ORS3.grmsAxis.all- ... 
config.ORS3.grmsAxis.low; 




ABC MPE UNNOTCHED ... 




[- 1 ,-I, conf ig. ABCUnnotched. ppod. X, conf ig. ABCUnnotched. force . X, . . . 

config.ABCUnnotched.control.X] = vibedataread(ABCUfilename, 'X+OdB. 
’); 

[- 1 ,-I, conf ig. ABCUnnotched. ppod. Y, conf ig. ABCUnnotched. force . Y, . . . 

config.ABCUnnotched.control.Y] = vibedataread(ABCUfilename, ’Y+OdB. 
'); 

[- 1 ,-I, conf ig. ABCUnnotched. ppod. Z, conf ig. ABCUnnotched. force . Z, . . . 

config.ABCUnnotched.control.Z] = vibedataread(ABCUfilename, 'Z-3dB. 
’) ; 

axis = [ ’xyz’ ]; 

for k = 1:length(axis) 

config.ABCUnnotched.ppod.Z.(axis(k)) = config.ABCUnnotched.ppod.Z.( .. 
axis(k))^2; 

end 

config.ABCUnnotched.force.Z = config.ABCUnnotched.force.Z^2; 

config.ABCUnnotched.control.Z = config.ABCUnnotched.control.Z^2; 
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%calculate Grms for each P-POD in every axis for every test axis 
config.ABCUnnotched.grms.X.all = grmscalc(20,2000^ ’p’ ^ . 

config.ABCUnnotched.ppod.X,4); 
config.ABCUnnotched.grms.X.low = grmscalc(20, cutoff p . 
config.ABCUnnotched.ppod.X,4); 

config.ABCUnnotched.grms.X.high = config.ABCUnnotched.grms.X.all- ... 
config.ABCUnnotched.grms.X.low; 

config.ABCUnnotched.grms.Y.all = grmscalc(20,2000, ’p’ , . . . 

config.ABCUnnotched.ppod.Y,4); 
config.ABCUnnotched.grms.Y.low = grmscalc(20,cutoff ,’ p . 
config.ABCUnnotched.ppod.Y,4); 

config.ABCUnnotched.grms.Y.high = config.ABCUnnotched.grms.Y.all- ... 
config.ABCUnnotched.grms.Y.low; 

config.ABCUnnotched.grms.Z.all = grmscalc(20,2000, 'p . 

config.ABCUnnotched.ppod.Z,4); 
config.ABCUnnotched.grms.Z.low = grmscalc(20,cutoff ,' p . 
config.ABCUnnotched.ppod.Z,4); 

config.ABCUnnotched.grms.Z.high = config.ABCUnnotched.grms.Z.all- ... 
config.ABCUnnotched.grms.Z.low; 

%calculate Grms for envelope of each P-POD 

ppodnum = 4 ; % # of P-PODs in data set - should be 8 unless it is .. 
unnotched ABC 

config.ABCUnnotched.max.PPOD = maxASD(config.ABCUnnotched.ppod, ... 

ppodnum); % max ASD values of each P-POD in each axis 

config.ABCUnnotched.grmsall.all = grmscalc(20,2000, 'p ',.. . 

config.ABCUnnotched.max.PPOD,4); 
config.ABCUnnotched.grmsall.low = grmscalc(20,cutoff, 'p ',.. . 
config.ABCUnnotched.max.PPOD,4); 

config.ABCUnnotched.grmsall.high = config.ABCUnnotched.grmsall.all- . 
config.ABCUnnotched.grmsall.low; 

%calculate Grms for envelope of each axis 

config.ABCUnnotched.max.Axis = maxAxis(config.ABCUnnotched.max.PPOD) 
%max ASD values of all P-PODs in each axis 
config.ABCUnnotched.grmsAxis.all = grmscalc(20,2000, ’d’ ,.. . 
config.ABCUnnotched.max.Axis,4); 
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config.ABCUnnotched.grmsAxis.low = grmscalc(20, cutoff^ 'd . 
config.ABCUnnotched.max.Axis,4); 

config.ABCUnnotched.grmsAxis.high = config.ABCUnnotched.grmsAxis.all- 
config.ABCUnnotched.grmsAxis.low; 

configl = config; 

save( 'testdata’configl’,’freq’,’-append’ ) % 


function grms = grmscalc(freq_start,freq_end,toggle,ppod,num) 

% function grms = grmscalc(freq_start,freq_end,toggle,ppod) 

%calculate Grms for a specified frequency range for all P-PODs 
%toggle is for # of P-PODs (’p’) or a single ASD 

global freq 

% freq_start = 20; 

% freq_end = 2000; 

index_start = find(freq==freq_start); 
index_end = find(freq==freq_end); 

axis = [ ’xyz’ ]; 

if toggle == ’p’ 
for k = 1:length(axis) 
for j = l:num; 

% for j = 1:size (ppod. (axis(k)),2); 

for i = index_start+l:1:index_end 

asdminusl = ppod. (axis(k)) (i-1,j); 

asd = ppod. (axis(k)) (i,j); 

db = lO^log(asd/asdminusl); 

oct = log(freq(i)/freq(i-1))/log(2) ; 

slope = db/oct; 

if 1 == index_start+l 

if slope == -3.010299957 

psdarea(j,k) = asdminusl^freq(i-1)^loglO(freq(1)/ 
freq(i-1)); 
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else 


psdarea(j^k) = (lO^log(2)^asd/(lO^log(2)+slope)) . 
■k (freq(i) - (freq(i-l) /freq (i) ) (slope/ (10^ log. . 
(2)))^freq(i-l)); 

end 

else 

if slope == -10^1og(2) 

psdarea(j^k) = psdareaminusl(j^k)+asdminusl^freq( 
1-1)^loglO(freq(1)/freq(1-1)); 

else 

psdarea(j,k) = psdareaminusl(j,k)+(lO^log(2)^asd . 
/(10^log(2)+slope))^(freq(1)-(freq(1-1)/freq( . 
1))^(slope/(10^1og(2)))^freq(i-l)); 

end 

end 

psdareaminusl(j,k)=psdarea(j,k); 

end 

% Grms(j^k) = (grms(end^j)); 

grms(j^k) = sqrt(psdarea(j^k)) ; 

end 

end 

else 

for k = 1:length(axis) 

for 1 = index_start+l:1:index_end 
asdminusl = ppod(i-l,k); 
asd = ppod(1,k); 
db = lO^log(asd/asdminusl); 
oct = log(freq(1)/freq(1-1))/log(2) ; 
slope = db/oct; 
if 1 == index_start+l 

if slope == -3.010299957 

psdarea = asdminusl^freq(1-1)^loglO(freq(1)/freq( 

i-1)); 

else 

psdarea = (lO^log(2)^asd/(lO^log(2)+slope))^ (freq 
(1)-(freq(1-1)/freq(1))^(slope/(10^1og(2))) . 
freq(1-1)); 

end 
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else 


if slope == -10^1og(2) 

psdarea = psdareaminusl+asdminusl^freq(1-1)^loglO 
(freq(i)/freq(i-l)); 

else 

psdarea = psdareaminusl+(lO^log(2)^asd/(lO^log(2) 
+slope))^(freq(1)-(freq(i-l)/freq(i))^(slope .. 
/(lO^log(2)))^freq(1-1)) ; 

end 

end 

psdareaminusl = psdarea; 

end 

% Grms(j^k) = (grms(end^j)); 

grms(k) = sqrt(psdarea); 

end 

end 


function maxPPOD = maxASD(config^ppodnum) 

%find max ASD values for each P-POD, including cross-axis responses 
testaxis = fieldnames(config); %pull test axis names 

axisnames = fieldnames(config.X); %only need one testaxis to grab ... 
axisnames 


for k = 1:numel(axisnames) % # of measurement axes 

for j = lippodnum % # of P-PODs 

for 1 = 1:numel(testaxis) % # of test axes 

dump(:,i,j) = (config. (testaxis{i}). (axisnames{k}) (:^ j)); 

end 

maxPPOD. (axisnames{k}) (: , j) = max(dump(:, : , j ) , [1^2); 
clear dump 
end 

end 


function maxAxis = maxAxis(config) 
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%find max ASD values for each axis^ including cross-axis responses 
fields = fieldnames(config); 
for k = 1:numel(fields) 

maxAxis(:,k) = max(config. (fields{k}), [],2)} 

end 


function [freq, a f^ppod,force^control] = vibedataread(filename^sheet) 

%W. Lan 
%4 DEC 2014 

%reads in 1 sheet of data from Excel, force summations, control ... 

channel, and outputs Grms for every P-POD in each 
%axis 

%read in 1 set of data (1 axis at 1 level) 

% [data,headers,all] = xlsread(’141118_NPSCuL+lsolators+ElatRing+8... 

P2M2s_WrongEL_Data.xlsx','MPE+6dB,X'); 

[data,headers,all] = xlsread(filename,sheet); 

% testaxis = input('Test axis (enter with single quotation marks, ... 
lower case): '); 

%pull out frequency 

freq = data(:,find(strcmp(headers(1,:),' X-Data [Hz]'))); 

%calculate frequency resolution 
Af = freq(2)-freq(1); 

%pull out control 

control = data ( : , find (-icellfun (Oisempty, strfind (headers (1, : ) , ' Control 
channel ')))); 

%find force summations (x=l, etc) 

f orce ( : , 1) = data ( : , find (-icellfun (Oisempty, strf ind (headers (1,:),'.. . 
Eorce Summation ' ) ) &-icellfun (Oisempty, strf ind (headers (1, : ) , 'X) ' ) ) ) . 

) ; 

f orce ( : , 2 ) = data ( : , find (-icellfun (@isempty, strf ind (headers (1,:),'.. . 
Eorce Summation ' ) ) &-icellfun (@isempty, strf ind (headers (1, : ) , ' Y) ' ) ) ) . 

); 
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force ( : ^ 3) = data ( : ^ find (-icellfun (@isempty^ strfind (headers (1 ^ ^ . 

Force Summation ' ) ) &-icellfun (@isempty, strf ind (headers (1, : ) , ' Z) ' ) ) 

); 


%find P-POD responses 

%P-POD DATA MUST BE ORDERED 1 THROUGH 8 

ppod. X ( : ^ : ) = data ( : ^ f ind (-icellfun (@isempty, strf ind (headers (1 ^ : ) ^ ' P 
POD ' ) ) &-icellfun (@isempty, strf ind (headers (1 ^ : ) , ' X ' ) ) ) ) ; 
ppod. y ( : , : ) = data ( : , find (-icellfun (@isempty, strf ind (headers (1, :) , ' P 
POD ' ) ) &-icellfun (@isempty, strf ind (headers (1, : ) , ' Y ' ) ) ) ) ; 
ppod. z ( : , : ) = data ( : , find (-icellfun (@isempty, strf ind (headers {1, :) , ' P 
POD ' ) ) &-icellfun (@isempty, strf ind (headers (1, : ) , ' Z ' ) ) ) ) ; 


clc; clear all; close all; 

%Post-processing script 

load testdata.mat 
% load isolator_trans_2OHz.mat 

global freq 

percentdifflv2a_ABCELVT.all = (configl.ABC.grmsall.all(l:6,:)-. .. 

config2a.ABC.grmsall.all)./configl.ABC.grmsall.all(l:6,:)^100; 
percentdifflv2a_ABCELVT.X = (configl.ABC.grms.X.all(1:6 . 

config2a.ABC.grms.X.all)./configl.ABC.grms.X.all(1:6,:)^100; 
percentdifflv2a_ABCELVT.Y = (configl.ABC.grms.Y.all(1:6,:)-... 

config2a.ABC.grms.Y.all)./configl.ABC.grms.Y.all(1:6,:)^100; 
percentdifflv2a_ABCELVT.Z = (configl.ABC.grms.Z.all(1:6,:)-.. . 
config2a.ABC.grms.Z.all)./configl.ABC.grms.Z.all(1:6,:)^100; 

"6 

%create overlay bar chart to show relative Grms of envelopes 
prettypictures(3)=figure (’units', 'normalized ', 'position', [0.0035 0. 

.0063 0.4930 0.9278 ],' Color ',[ 1 1 1]); 
highfreqpercent1 = configl.ABC.grmsAxis.high ./.. . 

configl.ABC.grmsAxis.all^lOO; 
widthl = 0.6; 
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% width2 = widthl/3; 

bar(highfreqpercent1^widthl^ 'FaceColor' ^ [1 1 1]) 
hold on 
grid on 

highfreqpercent2 = config2a.ABC.grmsAxis.high./.. . 

config2a.ABC.grmsAxis.all^100; 
bar(highfreqpercent2 ^ widthl) 

ylabel (' Percentage of G_{RMS} Above 500 Hz (%)',' FontSize 30) 
xlabel(’Test Axis', 'FontSize 30) 

set(gca, 'XTickLabel'FontSize' ,30) 
legendl=legend( 'Config 1' , 'Config 2' , 'Location' , 'SouthOutside 
Orientation','Horizontal') 
set(legendl, 'FontSize' ,30) 

% %loglog plots of data - max envelope by axis 
for j = 1:3 
% subplot (1,3,j) 

prettypictures(4 +j)=figure( 'units' , 'normalized', 'position', [0.5035 0. 

.0063 0.4930 0.9278 ],' Color ',[ 1 1 1]); 

% loglog(freq,configl.ABC.ppod.(test(i)).(axis(j))(:,position(j)),... 

freq,config2a.ABC.ppod. (test(i)) . (axis(j)) (:,position(j))) 
loglog(freq,configl.ABC.max.Axis(:,j),freq,config2a.ABC.max.Axis(: , j) 
) 

grid on 
hold on 

loglog (ones (2,1)^200, [lOE-7; lOEl] , ' g—' ) 

xlabel( 'Erequency (Hz) ', 'EontSize' ,30) 

dump = sprintf('ASD (g^2/Hz)'); 

ylabel(dump, 'EontSize' ,30) 

xlim([20 2000]) 

ylim([lOE-6 lOEl]) 

set(gca, 'EontSize ', 30) 

% str = sprintf('P-POD Response, %s Axis Envelope',test(j)); 

% suptitle (str) 

legendl=legend( 'Config l','Config 2', '200 Hz ',' Location ','.. . 

SouthOutside' , 'Orientation' , 'Horizontal' ) 
set(legendl, 'EontSize ', 30) 
end 
% end 
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% %loglog plots of data - max envelope by axis^ unnotched 
for j = 1:3 
% subplot (1,3,j) 

prettypictures(4 +j)=figure( ’units’ , ’normalized', 'position', [0.5035 0 
.0063 0.4930 0.9278 ],' Color ',[ 1 1 1]); 

% loglog(freq,configl.ABC.ppod. (test(i)) . (axis(j)) (:,position(j)), .. 

freq,config2a.ABC.ppod. (test(i)) . (axis(j)) (:,position(j))) 
loglog(freq,configl.ABCUnnotched.max.Axis(:,j),freq, ... 

config2a.ABCUnnotched.max.Axis(: , j)) 
grid on 

xlabel (’ Frequency (Hz) ', ’FontSize ', 30) 

dump = sprintf(’ASD (g'^2/Hz) ' ) ; 

ylabel(dump, 'FontSize’ ,30) 

set(gca, ’FontSize ’, 30) 

xlim([20 2000]) 

ylim([lOE-6 lOEl]) 

% str = sprintf('P-POD Response, %s Axis Envelope',test ( j)); 

% suptitle(str) 

legendl=legend( 'Config 1' , 'Config 2' , 'Location' , 'SouthOutside ',’... 

Orientation','Horizontal') 
set(legendl, 'EontSize’ ,30) 
end 
% end 


clc; clear all; close all; 

%Post-processing script 

load testdata.mat 
% load isolator_trans_2OHz.mat 

global freq 

% %loglog plots of data - max envelope by axis, unnotched 
for j = 1:3 
% subplot(1,3,j) 
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prettypictures(4 +j)=figure( ’units' , 'normalized', 'position', [0.5035 0 
.0063 0.4930 0.9278 ],' Color ',[ 1 1 1]); 

% loglog(freq,configl.ABC.ppod. (test(i)) . (axis(j)) (:,position(j)), .. 

freq,config2a.ABC.ppod. (test(i)) . (axis(j)) (:,position(j))) 
loglog(freq,config2d.ABC.max.Axis(:,j),freq, . . . 

config2d.ABCUnnotched.max.Axis(:,j)) 
grid on 

xlabel (' Frequency (Hz) ') 

dump = sprintf('ASD (g'^C/Hz) ' ) ; 

ylabel(dump) 

xlim([20 2000]) 

ylim([lOE-7 lOEO]) 

% str = sprintf('P-POD Response, %s Axis Envelope',test(j)); 

% suptitle (str) 

legend (' Config 3, ELVT' , 'Config 3, Unnotched ',' Location ','.. . 
SouthOutside' , 'Orientation','Horizontal','EontSize' ,20) 

end 
% end 


156 


APPENDIX G: 

MATLAB Code — Acoustic No-Test Rationale 

Calculation 


This appendix contains the MATLAB script for calculating the break-even surface areas 
to determine whether or not acoustic testing is required for P-PODs flying on an isolated 
NPSCuL. 


clc; clear all; close all; 

%Acoustic No-Test Rationale 
%W. Lan 
%9.30.2015 

fn = 100; %Hz - conservative assumption for panel mode 

PSD = 0.003; %G^2/Hz - lowest PSD value at 100 Hz from all isolated ... 
envelopes 

Q=10; %damping factor 

%use Miles Eqn for rms G's 
Grms = sqrt(pi/2^PSD^fn^Q); 

%response to acoustic noise - used OASPL from ABC User's Guide for 
%conservatism. Note: SPL at 100 Hz is 123.8 dB 
dB = 139.5; 

%Pressure Spectral Density 

Ps = 2.9E-9^2^10"'(dB/10) / (0.231^fn) ; 

%use Miles Eqn again nfor rms pressure 
Prms = sqrt(pi/2^Ps^fn^Q); 

%mass of payload - 6U or 2 3Us. Masses are low estimates to be ... 

conservative 
m_6U = 25; 
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m_3U = 35; 

%surface area calc 

SA_break_even_3U = Grms^m_3U/Prms 
SA_break_even_6U = Grms^m_6U/Prms 

SA_NPSCuL = 18.17^13 

break_even_ratio = SA_break_even_3U/m_3U %same as Grms/Prms 
% Following Eqns are not correct interpretation of "surface area" - 
not 

% surface area of component... should be surface area with Grms and . 
Prms 

% taken into account 

% break_even_ratio_3U = SA_NPSCuL/m_3U 
% break_even_ratio_6U = SA_NPSCuL/m_6U 
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APPENDIX H: 

Random Vibration Analysis Example Verification 


This appendix contains the MATLAB script and associated functions for calculating the 
frequencies, mode shapes, and response FRFs for a simply supported beam subject to a 
random vibration input at one end of the beam. The resulting response is compared with 
the results obtained using NASTRAN. The frequencies and mode shapes are consistent, 
but the response FRFs diverge after the first mode. This is most likely due to differences in 
effective modal masses between the two models. 


Table H.l: Frequency Comparison — MATLAB and NASTRAN FEMs 


Mode 

Frequency, Hz 

Percent Difference (%) 

MATLAB FEM 

NASTRAN FEM 

1 

30.9 

30.9 

0.267 

2 

124 

123 

0.941 

3 

279 

273 

2.04 

4 

495 

478 

3.53 

5 

774 

732 

5.36 
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Figure H.l: Mode Shapes of First Five Modes - Normalized Displacement 


160 









Figure H.2: Acceleration Response at L/4 and L/2 


clc; clear all; close all; 

% Modified from ME4522 PROJECT 3 
% W. LAN 
% 27 AUG 2015 

% global variables 

global E rho 1 A L M K ff x_vector 

% a. Definition of Nodal Coordinates 
elem=40; 

% elem=input('Number of elements: '); 
nodes=elem+l; 
q=nodes^2; 
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% c. Specification of material and geometric beam properties 

% material: aluminum 

E = le7; %lbf/in^2 

rho = 0.1; %lbm/in'^3 

rho=rho/386.4; %lbf-s^2/in^4 

% geometry 

L=10; %feet 

L=L^12; %inches 

x_vector=0:L/elem:L; 

W=5; %inches 
D=5; %inches 

% moment of inertia^ in^4 
I=W^D^3/12; 

% cross-sectional area, in^2 
A=W -k D ; 

% d. assembly of system stiffness and mass matrices 
k_global=zeros(q,q) ; 
m_global=zeros(q,q) ; 
for elem_index = 0:(elem-1) 
for i=l:4 

for j=l:4 

[ke me]=elem_ke(elem); 

k_global(2^elem_index+i,2^elem_index+j)=k_global(2^elem_index+i 
,2^elem_index+j)+ke(i,j); 

m_global(2^elem_index+i,2^elem_index+j)=m_global(2^elem_index+i 
,2^elem_index+j)+me(i,j); 
end 

end 

end 

%apply boundary conditions 

%For simply supported beam, DOFs 1 and n-1 are zero (left and right 
ends of beam) 

% new DOF 1 is rotation (old DOF 2) 

% keep only DOF 2 through q-2, q 

K=[k_global(2:q-2,2:q-2), k_global(2:q-2,q); ... 

k_global(q,2:q-2),k_global(q,q)]; 

M=[m_global(2:q-2,2:q-2) , m_global(2:q-2,q); . . . 
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m_global(q,2:q-2)^ m_global(q^ q) ]; 


%number of active DOFs 
ff=length(K); 

% HAPPENIN' NODES 
% node for L/4 
L4_node=nodes/4; 

L4_node=round(L4_node); 

% DOE for L/4 
L4_DOE=L4_node^2; 

% node for L/2 
L2_node=nodes/2-l; 

L2_node=round(L2_node) ; 

% DOE for L/2 
L2_DOE=L2_node^2; 

% g. Calculation of normal modes and natural frequencies 
format short g 
[v,d]=eig(K,M); 

% natural frequencies 
lambda=sqrt(nonzeros(d)) ; 
freq=lambda/(2^pi); %Hz 

omega = 20^2^pi:1^2^pi:2000^2^pi; 

% damped frequencies 
zheta=0.01; 

wd=lambda^sqrt(l-zheta^2); 
fd=wd/(2^pi); 

modal_M=v'^M^v; 
modal_K=v'^K^v; 
modal_C=diag(2^zheta^lambda); 

C=M^v^modal_C^v' 

% DOES for EREs 
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tf = [L4_D0F, L2_D0F]; 

for n = 1:length(omega) 

% INVERSE IMPEDANCE 

Z=[K-omega(n)^2^M+li^omega(n)^C] ; 

%transfer function 
H ( : , : , n) =inv (Z) ; 

end 

%normalize H 
for m = 1:ff 

for k = 1:ff 

h=abs(H(m,k,:)); 

H_norm(m,k, :)=h/(real(h (1,1,1))); 

end 

end 

figure 

subplot (2^1,1) 

loglog(omega,squeeze(abs(H(1,1,:))), '-b’ , ’LineWidth ’,1.5) 
subplot (2,1,2) 

loglog(omega,squeeze(H_norm(1,1, :)), '-k’ , ’LineWidth ’,1.5) 

%Random Input, Sx 

ABC_interface_MPE = [20^2^pi 40^2^pi 240^2^pi 2000^2^pi; 0.03 0.125 0. 
.125 0.003]; 

ABC_MPE_full = 10.^interpl(loglO(ABC_interface_MPE(1, :)) , loglO( . . . 
ABC_interface_MPE(2,:)),loglO(omega)); %ASD, g^2/Hz 

TE_drive = 2; 

%Apply input to left end 
Sx_vector=zeros(ff,length(omega)); 

% Acceleration at left end 
Sx=ABC_MPE_full; %g'^2/Hz 

Sx_vector (TE_drive, : ) =Sx_vector (TE_drive, : ) +Sx; 

% DOES for EREs 

tf = [L4_DOE, L2_DOF]; 
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for n = 1:length(omega) 

for m=l:length(tf) % number of tranfser functions 

Sxy(m,:,n) = abs(H_norm(tf(m),:,n)).^2^Sx_vector(:,n); 

end 

end 

[mode_shape_max, modeplot ] =modeshape (v) 

%load FEM results 

%FEM is constrained in translation and axial DOF at ends (simply ... 
supported) 

%FEM is constrained in 1 transverse direction for translation, axial 
%rotation and translation at all other nodes 
load Sxy_FEM_3.mat 
Sxy_FEM=Sxy_FEM_3; 

FREQ = omega/(2^pi); 

figure (’ Color ',[ 1 1 1]) 
subplot(2,1, 1) 

loglog(FREQ,Sxy(1, : ) , ’-b' , 'LineWidth ’,1.5) 
hold on 

loglog(FREQ,Sxy_FEM(:,1), '-k' , 'LineWidth ’,1.5) 
loglog(FREQ,ABC_MPE_full, '-g' , 'LineWidth' ,1.5) 
xlim([20 2000] ) 
grid on 

xlabel (' Frequency, Hz', 'FontSize' ,20) 

title (' Acceleration at Response Node, L/4' , 'FontSize' ,24) 
ylabel('ASD (g^2/Hz) ', 'FontSize ', 20) 
subplot (2,1,2) 

loglog(FREQ,Sxy(2,:), '-b' , 'LineWidth ',1.5) 
hold on 

loglog(FREQ,Sxy_EEM(:,2), '-k' , 'LineWidth ', 1.5) 
loglog(FREQ,ABC_MPE_full, '-g' , 'LineWidth ', 1.5) 
xlim([20 2000] ) 
grid on 

xlabel (' Frequency, Hz', 'FontSize' ,20) 

title (' Acceleration at Response Node, L/2' , 'FontSize' ,24) 
ylabel('ASD (g^2/Hz) ', 'FontSize ', 20) 
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legendl=legend( 'Matlab' , 'FEM'^ 'Random Input', 'Location', 'SouthOutside . 

' , 'Orientation' , 'Horizontal' ); 
set(legendl, 'FontSize' ,24) 


function [ke me] = elem_ke(elem) 

%output: k and m matrices for element 
%input: # of elements 
global E I rho A L 

l=L/elem; 

ke=E^I/1^3^ [12, 6^1, -12, 6^1;... 

6^1, 41 ^2, —6-^1, 2^1^2; . . . 

- 12 , - 6 ^ 1 , 12 , - 6 ^ 1 ; . . . 

6 ^ 1 , 2 ^ 1 ^ 2 , - 6 ^ 1 , 4 ^ 1 ^ 2 ]; 

me=rho^A^l/420^[156, 22^1, 54, -13^1;... 
22^1, 4^1^2, 13^1, -3^1'^2;... 

54, 13^1, 156, -22^1; ... 

-13^1, -3^1^2, -22^1, 4^1^2]; 


function [mode_shape_max,modeplot]=modeshape(v) 
global ff 

% first 5 mode shapes 

% separate displacement and rotation 
for 1=0: (ff/2-3) 

v_x(1+1,1:5)=v(2^1+2,1:5); 

% disp_x(1+1,1)=disp(2^1+1,1); 

% v_r(1 + 1,1:5)=v(2^1 + 2,1:5) ; 

end 

% v_x=[zeros(1,5);v_x;zeros(1,5) ]; 

% v_r=[zeros(1,5);v_r] ; 

modeplot=l; 

figure (' Color ',[ 1 1 1]); 
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% normalized mode shapes 
for i = 1:5 

mode_shape_max(:,!) = [v_x ( :,1) ./max(abs(v_x(:,!)))]; 

subplot (5^ 1, 1) 

plot(mode_shape_max(:^1)) 

grid on 

str = sprintf (' Mode %d',i); 
ylabel(str) 

% title (str) 

axis([l length(mode_shape_max) -1 1]) 
end 

xlabel( 'Active DOF') 
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